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I.  INTRODUCTION  AND  BACKGROUND 

The  transition  to  detonation  in  confined  granular  energetic  material  involves  a 
complicated  sequence  of  events  which  are  influenced  by  strength  of  confinement, 
particle  size,  permeability,  shock  sensitivity,  and  energy  release  rate  during 
combustion.  The  response  of  the  granular  bed  to  compressive  stress  waves  appears 
to  play  a  dominant  role  in  both  modes  of  deflagration-to-detonation  (DDT) 
transition  outlined  by  Bernecker,  Sandusky  k  Clairmont  [1]  and  Bernecker  [2J.  In 
Mode  I,  convective  combustion  near  the  ignitor  end  of  the  confined  column  creates 
a  rapidly  increasing  pressure  field  which  propagates  into  the  unburned  material  as 
a  compressive  wave  system.  This  stress  wave  system  can  compact  the  granular 
aggregate,  altering  the  propagation  rate  of  the  convective  flame.  If  the  pressure 
rise  from  combustion  is  sufficiently  rapid  to  drive  the  compaction  front  to  shock 
wave  strength,  then  unburned  material  may  be  ignited  by  mechanical  compression 
and  a  pathway  to  detonation  is  available.  In  Mode  II,  combustion  leads  to  a  mild 
rate  of  pressurization  of  the  bed,  although  compaction  of  the  aggregate  may  be 
extensive.  Suddenly  at  some  point  within  the  bed,  a  violent  reaction  drives  a  rapid 
pressure  build-up  (which  may  propagate  in  both  directions)  and  this  leads  to  shock 
wave  formation  and  provides  the  pathway  to  detonation.  In  these  scenarios,  the 
flow  of  gaseous  combustion  products  and  the  propagation  of  stress  waves  are 
intimately  coupled.  However,  an  initial  convective  flow  of  hot  gases  from  an 
ignitor  system  or  as  the  result  of  a  deflagration  is  not  a  prerequisite  for  transition 
to  detonation.  Impact  experiments  on  a  quiescent  granular  bed  show  that  high 
strain-rate  mechanical  compression  can  also  lead  to  a  transition  to  detonation. 

Special  laboratory  experiments  on  a  number  of  different  materials  have 
provided  many  important  observations  about  the  complex  transition  process.  The 
DDT  studies  of  Bernecker  k  Price  (3-5],  Price  k  Bernecker  [6],  Bernecker  et  aL 
[1,7,8],  the  combustion-driven  piston  experiments  of  Campbell  [9],  McAfee  k 
Campbell  [10],  MacAfee  et.  al.  [11],  and  the  shock-to-detonation  transition  (SDT)  or 
impact  studies  of  Sandusky  et.  al.  [12-15],  Green  et.  al.  [16],  and  Dick  [1748]  have  all 
contributed  to  our  understanding  of  the  behavior.  The  data  seem  to  suggest  that 
the  latter  stages  of  the  transition  are  essentially  independent  of  the  initial  stimulus, 
whether  from  impact  or  from  hot  gas  ignition.  Even  when  an  ignitor  system 
begins  a  convective  flame,  the  combustion  process  plays  the  role  of  a  piston  and 
drives  a  system  of  compression  waves  into  the  aggregate.  Thus,  successful 
transition  to  detonation  appears  to  require  formation  and  maintenance  of  a  strong 
compressive  wave  which  ignites  unburned  material  in  its  path.  The  present 
theoretical  studies  are  directed  toward  an  interpretation  of  this  wave,  before  a  full 
detonation  wave  has  been  established. 

Theoretical  predictions  of  the  transition  to  detonation  in  confined  granular 
energetic  material  have  been  underway  for  a  number  of  years.  Models  developed 
by  Baer  k  Nunziato  [19,20],  Weston  k  Lee  [21],  Aldis  et  al.  [22],  Price  et  al. 
[23-24],  Kim  [25],  Hsieh  k  Kim  [26],  Butler  k  Krier  [27],  Krier  k  Kezerle  [28], 
Akhatov  k  Vainshtein  [29,30],  and  Ermolaev  et  al.  [31]  originally  focused  on  DDT, 
but  most  now  have  the  ability  to  simulate  impact-generated  compaction  waves  also. 
These  analyses  employ  a  variety  of  constitutive  assumptions,  governing  equations, 
ignition  and  combustion  models,  and  numerical  solution  techniques.  These 
differences  can  lead  to  different  predictions  for  the  same  problem  as  illustrated,  for 
example,  in  Ref.  32.  However,  a  detailed  review  of  all  the  models  is  well  beyond 
the  scope  of  this  discussion.  Instead,  the  objective  of  this  report  is  to  summarize 
the  content  and  major  results  from  four  modeling  efforts  [33—36]  developed  recently 
by  the  present  author  in  an  attempt  to  understand  the  behavior  of  impact- 
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generated  compaction  waves.  Two  of  these  models  share  similar  features  with  the 
work  of  Baer  &  Nunziato  [19,20], 

The  ultimate  goal  for  any  model  is  predictive  capability,  although  few  ever 
achieve  it  A  less  spectacular  but  worthwhile  goal  is  the  simulation  of  an 
experiment,  where  ’’simulation"  may  involve  the  calibration  of  certain  model 
parameters  at  some  stage  in  the  comparison  or  calibration  to  data  from  similar 
runs.  Even  with  this  compromise,  the  model  simulation  can  sometimes  offer  an 
explanation  for  observed  behavior.  However,  before  attempting  any  simulation,  an 
absolute  necessity  is  the  ability  to  uniquely  specify  the  boundary  conditions  and/or 
source  terms  which  drive  the  experiment  If  the  simulation  problem  is  posed  such 
that  various  combinations  of  boundary  values  and  source  terms  could  be  chosen  as 
reasonable,  then  comparison  of  model  predictions  with  experimental  data  will  be 
meaningless.  A  major  problem  in  simulating  a  DDT  experiment  is  the  difficulty  in 
providing  a  unique  mathematical  description  of  the  ignitor  stimulus,  Le.,  specifying 

the  spatial  distribution  and  time-history  of  mass  flux  and  energy  content  of  the 

ignitor  system.  Impact  experiments  can  circumvent  this  problem  by  providing  a 
unique  time-dependent  boundary  condition  from  experimental  observations  of  the 
velocity  or  position  of  the  projectile  face  as  a  function  of  time.  The  present 

modeling  studies  were  designed  to  exploit  this  advantage. 

For  several  years  now,  the  Naval  Surface  Warfare  Center  has  been 
investigating  the  effects  of  impact  on  confined  granular  material  with  an  apparatus 
known  as  the  piston-driven-compaction  (PDC)  experiment  The  PDC  experiment 
[12-15],  shown  schematically  in  Fig.  1,  is  essentially  a  powder  gun  which  drives  a 
flat-faced  cylindrical  Lexan  projectile  into  a  confined  column  of  granular  material. 
Strength  of  confinement  can  be  varied  from  a  thick- wall  Lexan  tube  to  a  heavy- 
wall  steel  tube,  each  allowing  different  diagnostics.  Experimental  data  for  the 
location  of  the  projectile  as  a  function  of  time  after  impact  provide  a  unique 
boundary  condition  for  a  model.  Varying  the  initial  speed  of  projectile  impact 

produces  results  from  a  benign  compaction  wave  all  the  way  to  a  transition 
sequence  ending  with  high  order  detonation.  The  ability  to  predict  this  behavior 
would  enhance  the  credibility  of  any  theoretical  analysis. 

The  modeling  studies  summarized  here  are  focused  on  the  PDC  experiment 
The  objectives  of  this  effort  are  to  interpret  and  understand 

(a)  the  properties  of  the  initial  impact-generated  compaction  wave, 

(b)  the  compressive  ignition  behavior,  and 

(c)  the  acceleration  of  the  compaction  wave  to  a  strong  reactive  shock  wave. 

Each  of  the  four  models  to  be  discussed  will  attempt  to  address  one  or  more  of 
the  objectives  by  simulating  some  aspect  of  the  PDC  experiment,  often  within  the 
simplest  possible  framework.  Each  analysis  depends  in  some  way  on  a  constitutive 
model  of  the  granular  aggregate  (two-phase  mixture)  which  is  denoted  the 
"equilibrium  stress  state"  and  described  in  detail  in  Section  II.  Sections  HI— VI 
summarize  the  formulation  of  and  important  results  from  each  model: 

Section  III  -  Compaction  wave  behavior  with  no  reaction  or  combustion,  but 
porosity  adjustments  in  the  granular  aggregate  are  governed  by  a  rate-dependent 
mechanism. 
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Piston-Driven-Compaction  Experiment 
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Wall  Pressure  (Stress)  Time  History 


W//////////////////z/<^ 


Compacted 

Aggregate 


A 


Granular 

Material 


/ 

/ 

/ 

/ 

/ 

/ 


/ 


Compaction  Wave 


DIAGNOSTICS: 


Microwave 

Interferometry 


Lexan  Confinement:  Streak  Camera  Sc  Flash  Radiography 
Steel  Sc  Aluminum  Confinement:  Microwave  Interferometry, 
Wall  Transducers  &  Ionization  Probes 


Figure  1  -  Schematic  of  the  NSWC  piston-driven-compaction  (PDC)  experiment 
[12-15]  and  present  modeling  problem. 


Section  IV  -  Quasi-steady  compaction  wave  behavior  with  instantaneous 
adjustment  to  the  equilibrium  stress  state,  but  compaction/shock  wave  can  induce 
some  reaction  of  energetic  material  within  the  wave  front. 

Section  V  -  Simulation  of  collision  of  two  reactive  compaction/shock  waves 
(model  of  Section  IV)  along  with  predictions  for  the  resultant  flow  field  if  the 
collision  induces  an  increase  in  reactivity. 

Section  VI  -  Simulation  of  the  PDC  experiment  with  a  transient  reactive  shock 
wave  theory  which  accounts  for  wave-induced  reaction,  gas-phase  combustion  of  a 
reactive  intermediate  species,  and  ignition  and  transient  combustion  of  the  granular 
solid. 


It  may  be  appropriate  here  to  review  just  what  constitutes  a  "model"  in  the 
present  context  Although  the  four  analyses  are  different,  each  must  contain 

(a)  a  system  of  conservation  laws  or  balance  equations  based  on  a  consistent 
set  of  assumptions, 

(b)  a  set  of  constitutive  submodels  for  the  gas  phase,  the  solid  phase,  the 
two-phase  mixture,  all  phase  interaction  effects  (drag,  heat  transfer,  etc.)  which  are 
physical  processes  occurring  on  a  smaller  scale  than  will  be  resolved  by  the 
numerical  solution  to  conservation  laws,  prescriptions  for  ignition  and  combustion  of 
both  phases,  and 

(c)  a  numerical  integration  technique  which  will  minimize  artificial  distortion  of 
the  solution  being  sought 

Much  of  this  information  is  outlined  in  the  discussion  to  follow,  but  the  reader  is 
urged  to  consult  the  given  references  for  more  detail  on  each  model. 
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II.  DESCRIPTION  OF  THE  GRANULAR  AGGREGATE 


Predicting  the  behavior  of  compacted  granular  energetic  material  with  an 
analysis  which  resolves  micro-mechanics  at  the  grain  level  is  virtually  an  intractable 
problem.  Clearly,  some  compromise  is  required.  To  the  author’s  knowledge,  all 

current  theoretical  analyses  used  to  predict  the  transition  to  detonation  assume  that 
each  elemental  volume  contains  a  binary  mixture,  i.e.,  that  portion  of  the  volume 
not  occupied  by  the  solid  is  assumed  to  be  gas.  This  introduces  a  porosity 
variable  to  define  the  fraction  of  the  volume  occupied  by  the  solid,  and  its 

counterpart  for  the  gas: 


e  =  solid  volume  fraction  =  solid  volume  /  total  volume 
s  ' 


£ 


g 


1  - 


e  =  gas  porosity  =  gas  volume  /  total  volume 
s 


A  value  of  solid  volume  fraction,  however,  cannot  differentiate  among  particle  sizes 
or  shapes;  it  is  simply  a  state  property  associated  with  an  aggregate  of  ill-defined 
granular  material.  One  consequence  of  this  simplification  is  the  necessity  for  a 
phenomenological  sub-model  which  can  link  porosity  of  the  aggregate  to  its  stress 
state.  Of  course,  it  should  be  anticipated  that  the  value  of  certain  parameters  will 
be  determined  by  matching  the  model  response  to  experimental  data  for  the 
material  in  question.  In  a  similar  vein,  other  phenomenological  sub-models  will  be 
required  to  describe  the  drag  force  and  heat  transfer  between  phases.  These  must 
also  be  externally  calibrated  on  the  basis  of  experimental  data. 


A.  Implications  of  Constitutive  Assumptions 

The  particular  choice  of  a  sub-model  to  describe  the  dependence  of  the 
aggregate  stress  state  on  porosity  has  a  significant  influence  on  the  solution  to  the 
equation  system.  For  situations  of  interest  here,  the  constitutive  behavior  of  the 
pure  gas  phase  can  be  assumed  strain-rate  independent.  Similarly,  the  constitutive 
behavior  of  the  solid  at  theoretical  maximum  density  (TMD)  can  also  be  assumed 
rate  independent,  although  this  is  not  universally  true.  However,  since  compression 
of  an  aggregate  of  granular  solids  involves  simultaneous  sliding,  deformation,  and 
possibly  fracture  (depending  upon  the  type  of  material),  it  is  highly  likely  that 
rate-dependent  resistance  forces  will  be  encountered.  These  unknown  forces  are 
sometimes  modeled  with  a  global  parameter  called  compaction  viscosity  [see  Refs. 

19,  20,  33  and  Section  III]  in  a  rate-dependent  equation  which  relates  a  differential 
force  field  to  the  time  derivative  of  a  porosity  variable.  If  the  governing  equation 
system  for  the  mixture  neglects  diffusive  effects  such  as  Fourier  heat  conduction 
(not  interphase  heat  transfer)  and  viscous  stresses  (not  interphase  drag),  but 
includes  the  rate-dependent  description  of  porosity  change,  the  solution  will  predict 
compaction  waves  with  finite  thickness.  This  can  be  traced  to  the  dispersive  effect 
created  by  the  rate-dependent  porosity  equation.  For  most  compaction  waves,  the 
PDC  experiments  [12-15]  indicate  an  upper  bound  for  wave  thickness  which  is 

slightly  greater  than  6  mm.  If  an  analysis  employs  the  assumption  that 
compaction  viscosity  is  vanishingly  small,  then  porosity  instantaneously  adjusts  to  its 
new  equilibrium  value  and  the  rate-dependent  equation  is  replaced  by  an  algebraic 
constraint  equation.  Under  this  assumption,  the  solution  to  the  reduced  system  of 
governing  equations  predicts  that  compaction  waves  are  infinitesimally-thin  shock 
waves.  If  the  numerical  solution  predicts  otherwise,  the  finite  thickness  is  due  to 

artifical  diffusion  added  by  the  numerical  integration  technique.  The  modeling 

choice  between  rate-dependent  and  rate-independent  porosity  adjustment  creates 
much  the  same  situation  as  in  gas  dynamics:  the  inviscid  assumption  leads  to 
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shock  waves  as  discrete  discontinuities  while  retaining  diffusive  transport  leads  to 
shock  waves  with  finite  thickness.  The  model  discussed  in  Section  III  employs  a 
rate-dependent  porosity  adjustment  (similar  to  one  used  quite  successfully  by  Baer 
and  Nunziato  (19,20]),  while  the  models  discussed  in  Sections  IV- VI  employ  an 
instantaneous  (rate-independent)  porosity  adjustment  to  the  equilibrium  stress  state 
and  track  compaction  waves  as  moving  discontinuities. 


Another  important  consideration  is  the  treatment  of  the  solid  phase  material 
Compaction  of  a  granular  aggregate  can  involve  changes  in  density  of  the  solid 
particles  as  well  as  rearrangement  and  deformation  of  the  particles.  Quasi-static 
compaction  experiments  such  as  those  conducted  by  Elban  [37,38]  and  Campbell, 
Elban  &  Coyne  [37]  monitor  various  stresses  as  a  function  of  mixture  volume.  To 
deduce  the  associated  values  of  porosity  (or  solid  volume  fraction)  requires  a 
theoretical  model  of  the  mixture.  If  solid-phase  density  is  assumed  constant 
(incompressible),  then  porosity  values  are  determined  directly  by  the  observed 
volume  change.  For  materials  of  interest  in  this  study,  the  incompressibility 
assumption  can  easily  lead  to  predictions  of  mixture  density  which  have  exceeded 
the  value  of  TMD  at  atmospheric  pressure  (recall  that  TMD  is  a  function  of 
pressure  and  internal  energy).  Furthermore,  predictions  of  wave  propagation  speed 
in  the  mixture  will  be  in  error.  The  present  study  will  avoid  these  problems  by 
accounting  for  compressibility  of  the  solid-phase  material 


B.  Compressibility  of  the  TMD  Solid 


Compressibility  of  the  solid  phase  at  TMD  is  predicted  on  the  basis  of  an 
assumed  linear  path  in  the  shock  velocity  (V  )  /  particle  velocity  (u)  plane  as 
outlined  in  [40].  If  the  Hugoniot  of  the  TMD  solid  is  represented  as 


Vs  =  Sh 


+  A, 


sh 


u 


then  the  solid-phase  equation  of  state  can  be  written  as 


(1) 


P.  (/>s.e,)=^0fo[(eS-e>o)/(R+1)]+P.o[i+r0R/(R+i)] 

-(Psoash2/2)[2+roR/(R+1)]R/  (1+\hR)2  <2> 

where  R  =  ( PSq  /  P8)  -  1 


Prediction  of  pressure  versus  density  (Eq.  2  with  a  .  =  2.74  km/s,  A  .  =  2.65  and 
T  =  1.1)  for  an  isothermal  compression  of  TMD  nMX  is  shown  in  *Fig.  2  below 
afong  with  a  comparison  to  the  prediction  from  the  Helmholtz  Free  Energy 
function  used  by  Baer  &c  Nunziato  [19],  and  data  from  Olinger  et.  al.  [41].  The 
two  theories  give  almost  identical  results  which  compare  quite  favorably  to  Olinger's 
data. 
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theoretical  maximum  density  (TMD).  Data  for  0-HMX  are  from  Olinger  et  al.  [41). 
Solid  line  (P  )  from  Eq.  (2).  Dashed  line  (P.  )  is  from  theory  of  Baer  &  Nunziato 


C.  Equilibrium  Stress  State  of  the  Mixture 


The  solid  phase  material  exhibits  strength,  and  unlike  the  gas  phase,  can 
support  deviator  stress  components.  Anticipating  a  uniaxial  strain  calculation,  the 
two  principal  stresses  in  the  solid  phase  will  be  labelled, 

0g  =  axial  component  of  solid-phase  stress 

Or  =  radial  component  of  solid-phase  stress 

with  a  sign  convention  that  positive  values  are  in  compression.  Two  major 
assumptions  will  be  employed  here.  The  first  assumption  is  that  the  mixture 
stress  tensor  can  be  represented  as, 

a  =  €0  +  e_P_  =  axial  component  of  mixture  stress  tensor 

mx  ®  s  s  s 

0mr  =  eg0r  +  CgPg  =  radial  component  of  mixture  stress  tensor 

The  second  assumption  concerns  a  unique  relationship  between  0  and  0  which  is 
borne  out  by  the  experimental  compaction  data  [39,42,43]  for  several  materials  of 
interest  here.  This  is 


=  *  (Cs) 


where  £ 


mr 


mx 


experimental 


(4) 


which  is  illustrated  in  Fig.  3  below  for  the  double-base  ball  propellant  denoted 

TS-3659.  Note  that  under  conditions  of  the  quasi-static  compaction  experiment 

[37-39],  the  contribution  of  gas-phase  pressure  is  negligible  and  hence  the  definition 
of  {(eg)  closely  approximates  the  ratio  of  0_/0g.  The  values  of  solid  volume 
fraction  associated  with  the  data  points  in  Fig.  o  were  computed  from  experimental 
values  of  mixture  density  and  theoretical  estimates  of  solid-phase  density  based  on 
Eq.  (2).  The  solid  line  shown  in  Fig.  3  represents  the  function  assigned  to  {(c  ); 
note  that  the  aggregate  begins  deforming  as  a  fluid  medium  ({  =  1)  at  e  . 

When  the  mixture  deforms  in  uniaxial  strain,  the  spherical  stress  component, 

P  .  follows  directly  from 
m 

Pm  =  “mx  +  2t,mr>/3  =  UK  +  20  r>/3  +  egPg  <5> 

Invoking  the  relationship  in  Eq.  (4)  produces 

(®s  +  2V  !  3  =  °s  *e8)  where  *es>  =  l1  +  2««s>l  ’  3 

which  uniquely  defines  the  function  #(c  ).  If  the  spherical  component  of  the  solid 

phase,  P  ,  is  defined  as 
s 

ps  =  <®,  +  2V  /  3  =  os  Hes)  , 
then  the  mixture  pressure,  P  ,  can  be  written  as 

Pm  =  esPs“s’  es>  +  sgPg“g>  eg>  (7) 
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TMD.  Solid  line  is  function  {(c  ).  Experimental  data  are  values  of 
from  Run  #8  in  Ref.  42.  s 


Constitutive  theory  must  supply  prescriptions  for  P  ,  P  ,  and  a  method  to 
compute  €g.  The  phenomenological  model  of  porosity  change  i§  the  aggregate  is  a 
modified  version  of  the  concept  originally  presented  by  Carroll  &  Holt  [44].  Under 
quasi-static  conditions,  the  equilibrium  constraint  on  solid  volume  fraction  is  a  force 
balance  in  the  axial  direction  given  by 


(p.. 


£s)  -  pg  (pg  •  eg)  = 


Ps  (p»-  es) 

«(£») 


-pg  (Pg’  eg)  =M£s) 


(8) 


The  quasi-static  compaction  behavior  of  a  given  granular  mixture  is  built  into  the 
term,  0g,  which  represents  the  "strength"  of  the  aggregate  in  the  axial  direction 
which  must  be  exceeded  to  further  increase  the  solid  volume  fraction.  It  is 
represented  as 

ps  (£s)  =  Tl  [  1  -  fP'  +  B2  (f  P2  -  1  )  ]  ln  ( 1  /  £g)  (9) 

where  <*  s  (eg  /  £g0) 

and  Tj,  pj,  Bg,  and  pg  are  constants  to  be  determined  by  comparison  to 

experimental  data,  has  units  of  stress,  but  the  others  are  nondimensional. 

The  equilibrium  stress  state  is  determined  by  the  root  of  the  force  balance 
in  Eq.  (8).  Calibration  of  the  parameters  in  the  expression  for  0(eg)  [Eq.  9]  is  an 
iterative  process  to  produce  the  best  comparison  of  o  to  the  experimental  quasi¬ 
static  compaction  data  as  a  function  of  mixture  ciinsity.  The  results  for  an 

isothermal  compression  of  TS-3659  (Tj  =  2.6  Kpsi,  pj  =  5.5,  B,  =  32  and  p,  = 
0.08)  are  shown  in  Fig.  4.  The  dashed  line  represents1  theoretical1  values  of  mixture 
pressure,  P  ,  and  follows  from  Eq.  (7).  The  chain-dot  line  shows  the  variation  in 
TMD  predicted  by  Eq.  (2).  The  results  in  Fig.  5  for  the  radial  component  of 
mixture  stress  verify  the  utility  of  {(«  )  represented  with  the  solid  line  in  Fig.  3. 
Reference  34a  contains  an  appendix  which  lists  the  values  of  the  four  constants  in 
Eq.  (9)  for  a  number  of  different  propellants  and  two  inert  materials.  Note, 
however,  that  the  values  determined  in  Ref.  34  follow  from  the  early  analysis 
which  implicitly  assumes  that  the  mixture  stress  tensor  is  a  pressure  {%  =  1). 

Two  important  points  about  the  equilibrium  stress  state  deserve  emphasis. 
First,  the  two  phases  in  the  mixture  have  different  stress  states,  except  under 
conditions  of  a  fluidized  bed.  The  label  "equilibrium"  does  NOT  mean  that  both 
phases  are  forced  to  coexist  at  the  same  stress  level.  Instead,  equilibrium  implies 

the  final  relaxed  state  if  the  process  were  time-dependent  And  secondly,  the 

equilibrium  stress  state  is  calibrated  with  data  from  the  quasi-static  compaction 
experiment  [37-39]  under  conditions  where  the  contribution  from  gas-phase  pressure 
is  negligible.  However,  the  equilibrium  stress  state  follows  directly  from  the  root 
of  Eq.  (8)  which  is  indeed  sensitive  to  the  value  of  P  .  Combustion-generated  gas 
pressure  can  lead  to  an  equilibrium  stress  state  different  than  the  quasi-static 
state  for  equivalent  values  of  mixture  density.  Under  conditions  of  vigorous 
reaction  when  P  attempts  to  exceed  P  ,  tension  is  disallowed  by  setting  0g=O  and 
the  mixture  "unloads"  to  a  fluidized  bea 


9 


QUASI— STATIC  STRESS  STATE  *  60%  TMD  TS-3859  *  59%  Data:  Eban,Coyne, Campbell 
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Figure  5  -  Radial  component  of  mixture  stress  as  a  function  of  mixture  density 
for  an  isothermal  compression  of  6096  TMD  TS-3659  propellant.  Solid  line 
represents  a  _  from  Eq.  (3),  dashed  line  is  P  ,  and  chain-dot  line  is  TMD  from 
Eq.  (2).  Daia*  from  Run  #8  in  Sandusky  et.  al.  [42]. 


D.  Rate-dependent  Porosity  Change 


As  mentioned  in  Section  IIA,  compression  of  a  granular  aggregate  involves 
simultaneous  grain  sliding,  deformation  and  possibly  fracture  (if  the  material  is 
brittle),  and  hence  it  is  likely  that  strain-rate  dependent  resistance  forces  will  be 
encountered.  To  assess  the  importance  of  this  effect,  one  of  the  modeling  studies 
(Section  III)  provided  for  rate-dependent  porosity  change  of  the  aggregate.  In  a 
manner  similar  to  the  work  of  Baer  &  Nunziato  [19,20],  these  unknown  resistance 
forces  are  modeled  with  a  global  parameter  called  compaction  viscosity,  u  ,  which 
linearly  relates  the  time  derivative  of  solid  volume  fraction  along  the  solid-phase 
streamline  to  the  differential  force  field  describing  the  equilibrium  stress  state. 
Thus,  for  the  case  when  |  =  1, 

D(es)/Dts  =  (eg/ig  [Pg  -  Pg  -  0(*s)]  (10) 

The  influence  of  this  mechanism  on  compaction  wave  behavior  will  be  discussed  in 
Section  IIIA. 


E.  Gas-phase  Equation  of  State 


A  reactive  gas-phase  mixture  is  modeled  with  an  elementary  description 
which  envisions  two  gas-phase  species: 


and 


Yj  =  mass  fraction  of  reactive  intermediate  species  Rj  at  energy  level  Aegr° 
Yg  =  mass  fraction  of  final  combustion  product  Rg  at  energy  level  e  °. 


Mass  conservation  implies  Y,  +  Yp  =  1.  Section  VI  discusses  the  finite  rate  of 
conversion  from  R.  to  Rg  wnich  has  the  major  effect  of  simulating  the  release  of 
energy  as  intermediate  "reactants"  form  final  combustion  products.  The  gas-phase 
mixture  is  described  with  a  simple  constant  co-volume  equation  of  state  given  by, 


Pg0  /P«-b)  =  (R./J/)T,  (U) 

and  eg  =  cv  (lg  -  T°)  +  (y,  Aeffr°  +  Y2  eg°) 

where  y  (c  /c  ),  M  (molecular  weight),  and  b  (co-volume)  are  constants  whose 
values  are  Selected  as  reasonable  averages  from  BLAKE  Code  [45]  runs  for  each 
propellant  composition.  This  representation  only  approximates  the  real  gas  behavior 
of  the  mixture,  but  it  is  adequate  for  the  present  investigation  except  near 
conditions  of  detonation. 

F.  Phase  Interaction  Effects 


The  flow  of  gases  through  the  aggregate  creates  a  drag  force  per  unit 
volume,  f^,  which  is  modeled  with  the  Jones-Krier  correlation  [46], 
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(12) 


modified  by  the  tortuosity  factor  used  by  Gough  [47]  in  the  NOVA  interior 
ballistics  code.  Before  grain  burning  commences,  convective  heat  transfer  rates 
from  hot  gas  to  cold  solid  are  computed  from  the  convective  heat  transfer 
coefficient, 

hg  =  2kg(l.0  +  0.2Pr1/3Re2/3)  /dp  (13) 


Ignition  of  the  solid  is  assumed  to  occur  when  the  solid-phase  surface  temperature 
exceeds  an  arbitrary  constant  value  (e.g.,  473°K  for  TS-3659  ball  propellant).  The 
time-dependent  calculation  of  surface  temperature  is  done  with  Baer's  [48]  method 
which  approximates  the  solid-phase  thermal  profile  with  a  polynomial  and  creates 
an  ordinary  differential  equation  in  time.  Heat  transfer  and  combustion  enter  the 
two-phase  mixture  equations  as  "m"  and  "q"  which  are  per  unit  total  volume. 
Because  several  of  the  materials  of  interest  here  are  spherical  ball  propellants,  the 
present  analysis  has  extended  the  work  of  Jacobs  &  Sandusky  [49]  and  derived  an 
expression  for  S.  ,  the  exposed  surface  area  per  unit  total  volume  associated  with 
the  deformed  lattice, 
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"g(Rp/“p)~ 

_f(Rp/ap)_ 


(14) 


where  ap  =  distance  between  deformed  particle  centers 
Rp  s  radius  of  deformed  particle 

f(Rp/ap)  =  (2Iip/ap)3-no[(Rp/ap)-1/2]2[(4Rp/%)-'/2] 

6  (Rp  /  ap)  -  (2  Rp  /  ap)  {(2  Rp  /  ap)  “  nc  [(Rp  /  ap)  “  1  /  2]| 

Both  "m”  and  "q"  are  computed  on  the  basis  of  Sty. 

The  wall  boundary  friction  force  can  be  computed  from  data  provided  by 
the  quasi-static  compaction  experiment  [37,38]  as, 

friction  force  =  (applied  force- transmitted  force)/(y  x  diameter  x  length) 


Then, 

coefficient  of  boundary  friction  =  friction  force  /  o 

This  coefficient  for  TS-3659  ball  propellant  is  plotted  in  Fig.  5  of  Ref.  36,  and  is 

shown  to  increase  as  the  aggregate  is  compressed.  The  computations  for  TS-3659 

presented  here  have  assumed  the  worst  case  and  taken  the  maximum  value  of  0.20 

as  a  constant  Thus,  f  ,  =  0.20  a  ... 

’  w  mr 
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III.  COMPACTION  WAVE  MODEL  WITH  FINITE-RATE  POROSITY  CHANGE  [33] 


This  model  was  the  initial  effort.  The  objective  was  to  understand  the  early 
stages  of  the  PDC  experiment  by  simulating  the  propagation  of  an  inert 

impact-generated  compaction  wave  in  confined  granular  material.  Hence,  all 
combustion  and/or  reaction  events  are  deliberately  ignored.  As  in  all  the  models, 

deformation  of  the  aggregate  (two-phase  mixture  of  compressible  solid  grains  and 
gas)  occurs  in  uniaxial  strain.  In  this  model,  traction  forces  (friction)  between  the 
sliding  aggregate  and  the  lateral  confining  boundary  are  assumed  negligible.  Since 

the  quasi-static  compaction  experiment  was  not  yet  measuring  a  radial  stress 
component  when  this  model  was  under  development,  the  mixture  stress  tensor  is 

represented  with  the  spherical  component,  i.e.,  a  pressure  only  ({  =  1). 

The  central  focus  here  is  the  behavior  of  the  granular  aggregate.  As 

mentioned  in  Section  II,  compression  of  a  granular  mixture  involves  simulataneous 
grain  sliding,  local  grain  deformation  and  possibly  fracture  if  the  solid  material  is 
brittle.  It  is  quite  probable  that  these  events  will  give  rise  to  strain-rate 

dependent  resistance  forces.  An  assessment  of  this  influence  is  important  Because 
the  micro-mechanics  at  the  grain  level  present  an  intractable  computational 

problem,  the  combined  effect  of  these  unknown  forces  is  modeled  with  a  global 

parameter  called  compaction  viscosity,  u  .  Since  this  is  not  a  fundamental  material 
property,  its  magnitude  will  have  to  Le  determined  by  calibration  of  the  final 
solution.  Similar  to  the  formulation  used  by  Baer  &  Nunziato  [19—20],  compaction 
viscosity  linearly  relates  the  time  rate-of-change  of  solid  volume  fraction  along  the 
solid-phase  streamline  to  the  differential  force  field  describing  the  equilibrium  stress 
state.  This  relationship  is  given  by  Eq.  (10),  which  shows  that  the  equilibrium 
stress  state  is  recovered  as  the  volume  strain-rate  vanishes.  The  remaining  six 
balance  equations  (conservation  of  mass,  momentum  and  energy  for  each  phase)  are 

given  in  Ref.  33,  and  are  quite  similar  to  those  listed  in  Section  VL  The  total 

system  of  seven  equations  was  re-written  in  characteristic  form  (no  small  feat  in 
itself)  and  the  numerical  solution  was  obtained  by  the  method-of-characteristics. 

The  rate-dependent  porosity  adjustment  creates  a  dispersive  effect  on  the 
compaction  wave,  i.e.,  wave  thickness  will  tend  to  increase  with  distance  propagated. 
However,  nonlinear  convection  attempts  to  steepen  the  wave  profile,  counterbalancing 
the  dispersive  effect.  For  the  case  of  constant  velocity  impact  on  the  granular 
bed,  the  prediction  shows  a  compaction  wave  which  rapidly  approaches  a  constant 
thickness  and  steady  propagation  velocity,  in  accord  with  the  results  of  Baer  & 
Nunziato  [19].  The  original  calculations  in  Ref.  33  which  indicated  a  continuously 
spreading  wave  were  the  result  of  an  interpolation  error  in  the 

method-of-characteristics  solution.  Consider  an  example  case  of  constant  velocity 
(100  m/s)  piston  impact  on  a  quiescent  bed  of  7396  TMD  Class  D  HMX.  The 

(corrected)  prediction  for  compaction  wave  structure  is  shown  in  Fig.  6  as  a  plot 
of  mixture  pressure  versus  distance  at  20.5  us  after  impact.  In  this  example,  the 
compaction  viscosity,  u  ,  is  1000  poise  and  the  wave  thickness  is  approximately  5 
mm.  The  predicted  steady  propagation  speed  is  450  m/s,  while  Sandusky  observed 
438  m/s  [PDC-22,  Ref.  14];  the  difference  between  the  two  values  is  well  within  the 
variation  attributable  to  the  particle  size  distribution  of  Class  D  HMX. 
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Figure  6  -  Constant-velocity  piston  impact  at  100  m/s  on  7396  TMD  Class  D 
HMX.  Plot  shows  distribution  of  mixture  pressure  versus  distance  at  a  time  of 
20.5  (ta  after  impact  Compaction  viscosity  is  1000  poise.  Numerical  solution 
obtained  with  the  method-of-characteristics. 


There  are  several  important  observations  here.  The  first  is  that  the 
assumed  value  of  compaction  viscosity  influences  only  compaction  wave  thickness. 
As  intuition  might  suggest,  increasing  the  value  of  compaction  viscosity  creates  a 
"thicker"  compaction  wave.  In  the  PDC  experiment,  Sandusky  [13,14]  estimates  an 
upper  bound  on  wave  thickness  in  HMX  and  several  other  materials  as  the 
distance  between  tracer  wires,  6.35  mm,  which  were  seeded  into  the  bed.  The 
model  simulations  show  that  predicted  compaction  wave  thickness  will  exceed  this 
upper  bound  if  the  value  of  compaction  viscosity  is  much  larger  than  1000  poise. 
Futhermore,  the  simulation  in  Fig.  6  demonstrates  that  the  region  behind  the 
compaction  wave  (outside  the  finite  thickness)  remains  at  the  steady  value  dictated 
by  the  equilibrium  stress  state.  Secondly,  the  computed  value  of  steady 
propagation  speed  is  independent  of  the  value  of  compaction  viscosity.  Both 
observations  taken  together  imply: 

(a)  The  fundamental  properties  of  the  compaction  wave  (propagation  speed,  and 
downstream  values  of  pressure,  density,  and  porosity)  are  apparently  functions  of 
the  equilibrium  stress  state.  Hence,  a  complicated  time-dependent  numerical 
solution  should  not  be  required  to  predict  these  properties. 

(b)  At  least  for  quasi-steady  propagation  of  these  compaction  waves  (with  uq 
near  1000  poise),  the  influence  of  strain-rate  is  confined  to  the  region  within  the 
wave  thickness. 

An  additional  observation  is  the  "dual"  structure  of  the  wave,  similar  in  nature  to 
that  seen  in  an  "elastic/plastic"  material.  In  the  present  case,  impact  on  a 
medium  with  a  rate-dependent  equation  of  state  leads  to  a  rapidly  decaying 
transient  structure.  The  major  adjustment  forced  by  the  abrupt  change  in  particle 
velocity  becomes  the  slow  moving  "plastic"  wave  which  is  recognizable  in  Fig.  6  as 
the  compaction  front  near  10  mm.  However,  the  decaying  transient  structure  also 
leaves  behind  a  fast  moving  component  which  becomes  the  precursor  wave.  This 
precursor  is  visible  in  Fig.  6  near  40  mm  and  propagates  with  decreasing  (and 
almost  negligible)  amplitude  at  a  speed  slightly  less  than  2  mm/ us.  The  PDC 
experiment  [12-15]  found  no  evidence  for  a  precursor  wave  in  either  the  microwave 
interferometry  data  or  the  movement  of  the  tracer  wires.  If  this  wave  is  present, 
its  influence  is  apparently  small.  Further  discussion  on  the  evolution  of  the 
decaying  transient  structure  caused  by  the  initial  impact  can  be  found  in  Ref.  33, 
along  with  a  reference  to  the  early  work  of  B.  T.  Chu  as  reported  in  Ref.  50. 
One  final  result  from  the  model  predictions  for  the  transient  structure  is  that  the 
material  adjacent  to  the  projectile  face  undergoes  a  larger  value  of  volume  strain 
rate  than  does  material  which  traverses  the  fully-formed  compaction  wave.  The 
implication  is  that  if  the  compressive  ignition  phenomenon  is  rate  sensitive,  then 
material  adjacent  to  the  projectile  face  should  ignite  first  This  is  consistent  with 
most  of  Sandusky's  observations  from  the  PDC  experiment  [12-15]. 
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IV.  EQUILIBRIUM  SHOCK  WAVE  THEORY  [34] 


This  modeling  study  was  an  effort  to  exploit  an  important  result  from  the 
rate-dependent  model  (Section  III):  the  fundamental  properties  of  the  compaction 
wave  are  apparently  functions  of  the  equilibrium  stress  state.  This  suggests  the 
interesting  limiting  case  when  compaction  viscosity  is  equal  to  zero,  and  hence 
rate-dependent  porosity  change  is  replaced  by  an  instantaneous  adjustment  to  the 
new  equilibrium  stress  state.  In  this  limit,  a  finite-thickness  compaction  wave 
becomes  an  infinitesimally-thin  shock  wave,  and  the  details  of  the  wave  structure 
vanish  into  a  mathematical  discontinuity.  Although  a  description  of  the  wave 
structure  is  lost,  the  simplified  equation  system  is  easier  (and  faster)  to  solve. 
This  is  rather  analogous  to  the  inviscid  assumption  in  gas  dynamics  which 
approximates  a  shock  wave  as  a  discontinuity.  An  additional  objective  here  is  to 
deliberately  pose  the  simplest  possible  model  which  might  help  identify  controlling 
mechanisms  in  the  behavior  of  compaction  waves  in  granular  energetic  material. 
Such  an  attempt  is  always  haunted  with  the  prospect  of  neglecting  an  important 
part  of  the  physics  and  thus  fatally  distorting  the  predictions.  Although  the 
results  should  be  interpreted  with  due  caution,  simplicity  can  sometimes  provide 
insight. 


Some  interesting  questions  arise  if  we  relinquish  the  idea  that  compaction 
waves  are  inert  What  changes  in  behavior  would  be  expected  if  the  compaction 
wave  itself  caused  a  portion  of  the  solid  phase  material  to  react?  As  a 
compaction  wave  propagates  through  granular  energetic  material,  the  solid  volume 
fraction  is  forced  to  increase  abruptly  (over  a  time  interval  proportional  to  wave 
thickness  divided  by  speed).  If  the  strength  of  the  compaction  wave  increases,  the 
granular  aggregate  must  absorb  this  work  of  compaction  at  increasing  strain  rate. 
It  is  quite  improbable  that  this  energy  can  be  deposited  uniformly  throughout  the 
aggregate.  It  is  more  likely  that  the  energy  is  concentrated  in  thin  zones  near 
grain  boundaries  or  in  other  regions  of  material  weakness.  Hence,  above  some 
strain-rate  threshold,  the  temperature  of  a  small  amount  of  solid  material  may 

exceed  the  value  to  begin  reaction.  Even  if  the  reaction  fails  (e.g.,  quenching  by 
heat  loss  to  surrounding  cold  solid),  a  small  amount  of  solid  has  been  converted  to 
high  temperature  gaseous  combustion  products.  These  hot  gases  will  be  trapped  in 
the  available  pores  and  voids  of  the  collapsing  aggregate.  Thus  the  reactive 
compaction  wave  leaves  behind  a  compressed  aggregate  whose  available  porosity  is 
filled  with  combustion  gases  at  elevated  temperature.  Note  the  important 
distinction  that  the  combustion  gases  are  created  locally  by  the  wave  front,  as 

opposed  to  being  driven  into  the  aggregate  by  a  convective  flow.  This  concept 
provides  a  mechanism  which  can  transfer  the  compressive  work  done  at  high  strain 
rate  in  a  thin  compaction  wave  into  a  potential  source  of  ignition  distributed 

throughout  the  volume.  Furthermore,  the  gaseous  combustion  products  trapped  in 

the  porosity  of  the  aggregate  can  also  effect  the  mechanical  "stiffness"  of  the 
compressed  aggregate,  and  hence  alter  the  rate  of  propagation  of  a  compaction 
wave.  Recently,  Dick  [18]  re-examined  data  from  his  earlier  SDT  experiments  on 

6596  TMD  HMX  and,  independent  of  this  work,  came  to  a  similar  conclusion  about 
the  existence  of  wave-induced  reaction. 

Compaction  waves  in  this  model  are  treated  as  reactive  shock  waves.  The 

mechanism  responsible  for  wave-induced  reaction  of  some  portion  of  the  solid 
material  is  undoubtedly  related  to  mechanical  work  done  at  high  strain  rate.  A 
description  of  this  process  in  a  granular  aggregate  involves  extremely  complex 

physics  and  is  presently  unavailable.  Of  course,  a  shock  wave  model  (zero 

thickness,  and  infinite  strain  rate)  could  not  predict  the  strain-rate  history  imposed 
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on  material  which  passes  through  the  compaction  wave.  To  circumvent  this 
problem,  the  analysis  assumes  that  a  prescribed  fraction  of  solid  material  reacts 
within  the  shock  front  to  form  gaseous  combustion  products  immediately 
downstream.  Outside  the  shock  front,  reaction  ceases.  In  reality,  reaction  probably 
does  fail  when  the  system  is  near  the  limiting  condition.  For  those  systems  which 
exceed  this  limit,  it  is  more  likely  that  reaction  "falters"  during  an  induction  period 
and  then  resumes. 

The  equilibrium  shock  wave  model  [34]  of  compaction  waves  intentionally 
neglects  all  "loss"  mechanisms  as  well  as  rate-dependent  material  processes. 
Resistance  to  porosity  change  is  assumed  vanishingly  small,  i.e.,  the  compressed 
aggregate  instantaneously  adjusts  to  its  equilibrium  stress  state  [see  Section  IIC]. 
Deviator  components  of  the  mixture  stress  tensor  are  ignored  and  the  equilibrium 
stress  state  is  based  on  a  mixture  pressure.  It  is  assumed  that  both  phases  travel 
at  the  same  velocity  in  the  compressed  aggregate,  which  eliminates  gas/particle 
drag.  Inter-phase  heat  transfer  is  also  neglected,  although  the  two  phases  may 
exist  at  different  temperatures.  The  heat  transfer  problem  is  addressed  separately 
in  Ref.  34,  and  this  artificial  decoupling  is  removed  in  Ref.  36.  Friction  between 
the  compressed  aggregate  and  the  lateral  confining  boundary  is  ignored. 
Combustion  occurs  only  within  the  wave  front  itself.  A  specified  fraction,  rj,  of 
the  granular  solid  material  can  react  and  form  gas-phase  combustion  products 
(with  a  given  energy  release)  during  passage  of  the  shock  wave,  i.e.,  shock-induced 
reaction.  Although  the  eventual  goal  is  to  be  able  to  predict  the  value  of  the 
important  parameter  r),  at  present  its  value  must  be  specified.  Thus,  in  certain 
cases,  this  analysis  will  be  used  to  interpret  experimental  data  instead  of  predicting 
it. 


The  model  problem  illustrated  in  Fig.  7  is  used  to  simulate  the  PDC 
experiment.  The  computational  problem  involves  two  shock  waves  and  a  contact 
discontinuity,  each  of  which  is  treated  as  a  moving  internal  boundary.  One  shock 
wave  (Vgj)  travels  back  up  the  Lexan  rod  (projectile)  to  arrest  its  initial  motion, 
and  the  other  (V  )  represents  the  compaction  wave  in  the  granular  material.  The 
contact  discontinuity  is  the  interface  boundary  between  the  compressed  Lexan  rod 
and  the  compacted  aggregate.  With  the  assumptions  outlined  above,  constant 
velocity  impact  leads  to  steady  propagation  of  all  three  waves;  there  are  no 
transients.  This  provides  a  substantial  simplification  in  the  solution  procedure 
which  is  discussed  in  Ref.  34  along  with  the  equation  system. 

It  is  important  to  note  how  the  assumptions  used  in  the  theory  enter  the 
comparisons.  All  "loss"  processes  have  been  neglected,  and  the  compressed 
aggregate  behind  the  compaction  shock  wave  is  constrained  to  the  equilibrium  stress 
state.  When  predictions  of  compaction  wave  speed  are  plotted  versus  projectile 
speed,  this  incorporates  the  assumption  of  an  ideal  momentum  transfer  at  projectile 
impact,  i.e.,  no  losses  at  impact  On  this  plot,  it  is  common  that  the  experimental 
data  lie  on  or  below  the  inert  shock  wave  prediction.  When  the  predictions  are 
plotted  versus  interface  speed  (called  particle  velocity  in  Sandusky's  work),  all  losses 
associated  with  projectile  impact  have  been  removed  from  the  comparison.  If 
experimental  data  lie  above  the  inert  shock  wave  prediction  in  this  plot,  some 
mechanism  (such  as  heat  release  during  reaction)  is  driving  the  wave,  or  dynamic 
compaction  of  the  aggregate  has  produced  a  stress  state  greater  than  the 
equilibrium  value. 
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Figure  7  -  Schematic  of  computational  problem  used  to  simulate  the  NSWC 
piston-driven-compaction  experiment  [12-15].  Compaction  wave  in  two-phase  mixture 
is  reactive  equilibrium  shock  wave. 
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Predictions  from  the  equilibrium  shock  wave  model  have  been  compared  to 
experimental  data  from  the  PDC  experiment  for  many  different  granular  materials 
(inert  and  energetic)  with  widely  varying  properties.  The  reader  is  invited  to 
consult  Ref.  34  for  details  on  a  number  of  these  comparisons.  Three  example 
cases  will  be  discussed  here. 

A.  #20  Sieve  Cut  from  Class  D  HMX 

Cyclotetramethylenetetranitramine  (HMX)  is  a  brittle  crystalline  energetic 
substance  which  is  a  major  ingredient  in  secondary  explosives  and  high  energy 
propellants.  Class  D  HMX  is  characterized  by  a  broad  particle  size  distribution, 
and  thus  successive  samples  taken  from  this  material  may  exhibit  different 
mechanical  properties.  To  eliminate  this  uncertainty,  Elban  &  Chiarito  [38]  prepared 
a  sieve  cut  from  Class  D  material  which  has  an  average  particle  size  of  925  /zm 
within  a  narrow  distribution.  Predictions  for  compaction  wave  speed  in  this 
material  at  initial  density  of  6096  TMD  are  shown  in  Fig.  8  versus  interface  speed 
(particle  velocity).  Figure  8  also  includes  data  from  two  PDC  experiments  [14] 
conducted  in  the  intermediate  confinement  provided  by  the  thin-wall  aluminum  tube 
which  eliminates  the  uncertainty  associated  with  radial  deformation  of  Lexan  tubes. 
The  good  comparison  between  theory  and  experiment  suggests  that  compaction 
waves  in  #20  sieve-cut  HMX  behave  as  inert  waves,  at  least  up  to  a  particle 
velocity  of  125  m/s.  Note  that  fracture  of  HMX  grains  during  compression  of  the 
aggregate  is  implicity  included  in  the  equilibrium  stress  state,  but  otherwise  has 
been  ignored.  Of  course  as  these  waves  increase  in  strength,  they  will  not  remain 
inert.  But  the  simple  theory  has  provided  an  important  base  line  to  assess  the 
behavior  of  stronger  waves. 

B.  WC-231  Propellant 

A  second  example  is  provided  by  WC-231  which  is  a  double-base  (25%  NG 
in  NC)  propellant  whose  grains  are  shaped  like  "fat  pancakes"  (rolled  spheres)  with 
an  average  diameter  of  approximately  790  microns.  Predictions  for  compaction 

wave  speed  in  50%  TMD  WC-231  versus  interface  speed  are  compared  in  Fig.  9  to 
data  from  several  of  NSWCs  PDC  and  CGC  (cold  gas  compaction)  experiments 
[13,14].  The  comparison  shows  that  the  experimental  compaction  wave  speeds  are 
greater  than  the  model  predictions  for  an  inert  wave.  Since  loss  mechanisms  have 
been  ignored,  something  has  to  be  "driving"  these  waves.  In  principle,  rate- 

dependent  resistance  to  dynamic  compaction  of  the  granular  aggregate  could  create 
an  increased  stress  state  which  would  support  the  higher  speed;  however,  rate- 
dependent  compaction  models  [19,20,  and  Section  III]  suggest  that  the  magnitude  of 
relaxation  time  required  to  maintain  this  condition  for  the  duration  of  the 

experiment  would  produce  excessively  thick  compaction  waves  (not  indicated  by  the 
experiment).  The  other  possibility  is  wave-induced  reaction  involving  a  small 
fraction  of  solid  material  (also  postulated  by  Dick  [18]).  Since  the  aggregate  is 
simultaneously  being  compressed,  the  newly  created  combustion  gases  released  into 
the  available  porosity  of  the  aggregate  are  trapped  and  will  exert  a  pressure.  This 
gas  pressure  increases  the  "stiffness"  of  the  two-phase  aggregate,  which  in  turn 
supports  an  increased  compaction  wave  speed.  Computations  show  that 

wave-induced  reaction  (which  releases  the  total  heat  of  reaction  of  the  solid)  of  a 
few  tenths  of  one  percent  of  the  solid  phase  creates  enough  gas  pressure  to  stiffen 
the  aggregate  and  enable  it  to  support  the  compaction  wave  speeds  measured  in 
the  experiments.  Sandusky's  [14]  comments  concerning  the  extent  of  reaction 
eventually  observed  in  two  of  the  experiments  are  noted  in  Fig.  9;  they  support 
(but  do  not  prove)  the  concept  of  a  reactive  compaction  wave. 

(text  continues  on  page  24) 
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Figure  8  -  Predicted  compaction  wave  speed  versus  interface  speed  for  6096  TMD 
#20  sieve  cut  HMX.  Solid  line  represents  inert  compaction  wave.  PDC  data  from 
Sandusky  [13]. 
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Figure  10  -  Predicted  compaction  wave  speed  vs  interface  speed  (particle  velocity) 
for  6096  TMD  TS-3659  ball  propellant.  Data  from  the  NSWC  PDC  experiment  [15]. 
Solid  line  is  inert  compaction  wave,  dashed  and  chain-dot  lines  involve 
wave-induced  reaction. 


C.  TS-3659  Ball  Propellant 


A  third  example  is  the  granular  propellant  used  as  the  basis  for  three  out 
of  four  modeling  problems  in  the  1988  JANNAF  Propulsion  Systems  Hazards 
Subcommittee  Modeling  Workshop  [32].  TS-3659  is  a  double-base  (21.5%  NG  in  NC) 

ball  propellant  with  nearly  spherical  grains  of  434  urn.  average  diameter. 
Predictions  for  compaction  wave  speed  are  compared  in  Fig.  10  to  NSWC  data 
from  the  three  PDC  experiments  which  backed-up  the  modeling  workshop  problems. 
As  in  the  case  of  WC-231,  the  PDC  data  lie  above  the  prediction  for  an  inert 
compaction  wave  suggesting  that  wave-induced  reaction  was  involved.  However,  the 
trend  in  WC-231  that  increasing  particle  velocities  (stronger  waves)  are  accompanied 
by  a  monotonic  increase  in  the  amount  of  wave-induced  reaction  is  reversed  in 
TS-3659.  This  reversed  trend  also  applies  to  WC-140  single-base  propellant  and 
ABL-2523  casting  powder.  The  comparison  in  Fig.  9  for  TS-3659  suggests  that  as 
the  waves  become  stronger,  a  smaller  percent  of  the  solid-phase  material  entering 
the  wave  front  is  involved  in  reaction.  This-  does  not  necessarily  mean  that  very 
strong  compaction  waves  will  behave  as  inert  shock  waves;  experimental  evidence  is 
quite  the  opposite.  One  possible  explanation  is  that  the  high  strain  rate  in  the 
stronger  waves  is  momentarily  creating  high  temperature  material  at  grain 
boundaries  which  becomes  viscous  or  liquid-like.  This  viscous  layer  permits  grain 
sliding  and  allows  the  aggregate  to  collapse  toward  the  equilibrium  stress  state  for 
inert  material,  without  a  proportional  increase  in  reactivity.  The  existence  of  these 
opposite  trends  (apparent  reactivity  increases  or  decreases  with  wave  strength) 
serves  to  underscore  how  little  is  known  about  the  micro-mechanics  at  the  grain 
level,  and  exactly  what  mechanisms  are  responsible  for  compressive  ignition  of 
material.  Clearly,  this  problem  will  require  further  study. 


V.  COLLISION  OF  REACTIVE  COMPACTION/SHOCK  WAVES  [35] 

It  is  common  practice  to  display  data  from  impact  experiments  (e.g.,  PDC 
[12-15]  at  NSWC),  DDT  experiments  (e.g.,  Bernecker  [8]  at  NSWC)  and  the 

combustion-driven  piston  experiment  (e.g.,  McAfee  [10,11]  at  LANL)  in 

"Distance-Time"  plots.  In  many  cases,  these  X-T  plots  appear  to  show  abrupt 
changes  in  propagation  speed  as  the  compaction  wave  accelerates  to  a  reactive 
shock  wave  (speeds  in  the  range  of  1-3  km/s)  and  then  to  a  detonation  wave. 
Furthermore,  Bernecker's  DDT  experiment  [8]  with  TS-3659  double-base  ball 
propellant  shows  the  presence  of  a  rearward  propagating  compressive  wave,  i.e.,  a 
wave  which  propagates  back  into  material  which  is  already  burning.  The  origin  of 
this  wave  seems  to  be  in  the  vicinity  of  an  increase  in  speed  of  a  strong 
compaction  wave. 

The  PDC  experiments  have  led  to  the  observation  [14]  that  the  time  to  the 
onset  of  compressive  reaction  correlates  as  "r  At  =  constant",  where  r  is  a 
measure  of  the  stress  state  in  the  solid  phase  and  At  is  the  time  interval  after 
impact.  Sandusky  notes  that  the  first  appearance  of  "flame"  usually  occurs  in 
material  adjacent  to  the  projectile  face.  Vigorous  combustion  in  this  region  will 
generate  and  support  a  stress  wave  system  which  must  propagate  into  material 
already  compressed  by  the  initial  compaction  wave.  The  objective  of  this  effort 
was  to  determine  if  collision  of  this  wave  system  with  the  leading  compaction 
front  might  produce  an  altered  wave  pattern  responsible  for  the  abrupt  change  in 
speed  during  acceleration. 
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A.  Model  Development 


Simulation  of  the  collision  between  two  compaction/shock  waves  of  the  same 
family  (the  weaker  wave  is  overtaken  from  behind)  employs  the  reactive  equilibrium 
shock  wave  model  discussed  in  Section  IV.  The  collision  model  follows  from  a 
simple  idea.  To  begin,  consider  a  piston  in  contact  with  a  quiescent  inert  rate- 
independent  medium  as  sketched  in  Fig.  11.  At  time  zero,  the  piston  assumes  a 
constant  velocity.  A  plot  of  this  event  in  the  "Distance-Time"  plane  shows  a 
conbtant  velocity  shock  wave  running  ahead  of  the  piston  path. 

Note  the  direct  analogy  of  this  wave  diagram  with  the  flow  field  produced 
by  steady  two-dimensional  supersonic  flow  over  a  wedge  (i.e.,  the  piston  path). 

Focusing  on  the  analogy  with  supersonic  flow  for  a  moment,  imagine  that  the 

single  wedge  is  replaced  with  a  double  wedge,  ie.,  the  body  has  a  second  abrupt 
change  in  slope  (see  shaded  region  in  Fig.  12).  Of  course,  this  new  slope  change 

will  create  a  second  shock  wave  which  eventually  must  collide  with  the  leading 
shock  wave.  As  illustrated  in  Fig.  12  for  a  nonreacting  medium,  the  collision 
produces  a  combined  shock  wave  whose  strength  is  intermediate  between  the 

strengths  of  the  two  waves  attached  to  the  body.  A  centered  simple  wave  system 
(expansion  fan)  anchored  at  the  collision  point  will  attempt  to  equilibrate  this 
pressure  difference. 

Of  course  steady  supersonic  flow  over  a  double  wedge  is  not  the  problem  of 
interest.  However,  application  of  this  wave  diagram  to  one-dimensional  wave 
propagation  requires  an  interpretation  of  the  second  increase  in  piston  velocity. 
Once  the  projectile  impacts  the  the  granular  bed  in  the  PDC  experiment,  it  does 
not  suddenly  accelerate  to  an  increased  speed  after  a  certain  time  interval.  On 
the  contrary,  its  velocity  may  begin  to  decrease.  However,  the  onset  of  combustion 
of  material  adjacent  to  the  piston  face  might  generate  a  sudden  increase  in  stress 
level  which  would  have  the  same  effect  as  an  increase  in  piston  velocity.  Thus, 
the  present  model  will  employ  an  increase  in  piston  velocity  to  simulate  the  effect 
of  vigorous  combustion  which  begins  at  the  piston  face,  after  a  suitable  delay  time, 
and  then  propagates  away  from  the  piston  into  the  compressed  aggregate.  Note 
that  the  combustion  process  which  "begins"  at  the  piston  face  is  the  onset  of 
substantial  burning  in  the  aggregate  since,  in  the  model  of  Section  IV,  a  small 
amount  of  the  solid  material  could  have  already  reacted  within  the  initial 

compaction  wave. 

The  wave  diagram  for  a  double  impact  shown  in  Fig.  12  applies  to  an  inert 
medium.  If,  however,  the  medium  is  reactive  and  shock  waves  are  capable  of 

inducing  energy  release,  the  wave  pattern  can  easily  be  altered.  In  particular,  if 
the  combined  shock  wave  induces  a  nominal  increase  in  reactivity  compared  to  the 
initial  shock  wave,  this  will  raise  the  downstream  mixture  pressure  (or  stress  level) 
enough  to  wipe  out  any  expansion  wave.  Any  further  increase  in  reactivity  will 
raise  the  mixture  pressure  such  that  a  rearward  propagating  shock  wave  will 
appear  in  the  location  of  the  original  expansion  wave  system,  as  illustrated  in  Fig. 
13.  The  dividing  streamline  (chain-dot  line)  which  passes  through  the  collision  point 
becomes  a  contact  discontinuity  between  the  regions  downstream  of  shock  waves  (3) 
and  (4)  (see  Fig.  13).  A  contact  discontinuity  separates  two  regions  with  different 
density  and  temperature,  but  particle  velocity  and  mixture  pressure  must  be  equal. 
The  physical  meaning  is  the  following.  All  the  material  lying  above  the  contact 

discontinuity  (closest  to  shock  wave  4)  has  been  influenced  only  by  shock  wave  (4); 
all  the  material  lying  below  has  been  subjected  to  shock  waves  (1),  (2)  and  (3). 

This  contact  discontinuity  plays  an  important  role  in  the  examples  to  follow. 

(text  continues  on  page  29) 
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Figure  11  -  Distance-time  plot  of  constant-velocity-piston  impact  on  a  quiescent 
inert  rate-independent  medium.  Note  analogy  with  steady  -  two-dimensional 
supersonic  flow  over  a  wedge. 
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Distance 


Figure  12  -  Distance-time  plot  of  double  impact  of  piston  (two  abrupt  changes  to 
constant  velocity)  on  a  quiescent  inert  rate-independent  medium.  Note  analogy 
with  steady  supersonic  flow  over  a  double  wedge. 
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Figure  13  -  Possible  distance-time  plot  of  double  piston  impact  on  quiescent 
rate-independent  reactive  medium. 
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The  problem  to  be  solved  is  that  posed  in  Fig.  13.  Both  piston  velocities 
are  assumed  constant.  Since  all  rate  processes  have  been  removed  from  the 
equation  system,  each  shock  wave  is  quasi-steady  and  travels  at  a  constant 
velocity.  Thus  the  solution  can  be  found  from  a  system  of  coupled  "jump 
condition"  equations  which  are  listed  in  Ref.  35  along  with  a  description  of  the 
solution  procedure.  These  equations  are  based  on  the  same  set  of  assumptions  as 
the  equilibrium  shock  wave  theory  of  Section  IV  and  will  not  be  repeated  here. 
Note  again  that  the  analysis  cannot  predict  the  amount  of  reactivity,  i),  in  the 
wave-induced  reaction;  this  value  must  be  deduced  from  comparison  to  experimental 
wave  patterns.  In  some  examples,  it  will  be  assumed  that  the  shock-induced 
reaction  releases  the  total  heat  of  reaction  of  the  energetic  solid.  However,  in 
single-  and  double-base  propellants,  the  available  evidence  suggests  that  only  a 
fraction  of  the  total  heat  of  combustion  is  released  in  this  reaction  (see  discussion 
below). 

B.  Results 


This  analysis  will  be  used  to  interpret  data  from  Sandusky's  [14,15]  PDC 
experiment,  a  combustion-driven  piston  experiment  of  McAfee  [10,11],  and  a  DDT 
experiment  reported  by  Bernecker  [8J.  Granular  materials  include  two  double-base 
ball  propellants  (WC-231  and  TS-3659)  as  well  as  Class  D  HMX.  Values  of  the 
modeling  parameters  used  to  define  the  quasi-static  stress  state,  total  energy 
release,  properties  of  the  gas-phase  combustion  products,  etc.,  for  these  materials 
can  be  found  in  Refs.  34a  and  35b. 

PDC-80  /  60.296  TMD  TS-3659  Ball  Propellant 


TS-3659  is  a  double  base  ball  propellant  (nominally  21.6%  NG  in  NC)  with 
essentially  spherical  grains  (434  Mm  diameter).  The  PDC  experiment  [15]  denoted 
PDC-80  involved  a  160  m/s  impact  on  602%  TMD  propellant  confined  in  the 
aluminum  tube,  101.7  mm  long.  Although  later  runs  in  this  series  made  use  of 
microwave  interferometry  to  track  motion  of  the  strong  compaction  front,  PDC-80 
relied  on  two  wall-mounted  pressure  transducers  located  at  38.1  and  762  mm  from 
the  impacted  end.  The  time  at  which  the  compaction  front  passes  the  gage 
location  is  determined  as  the  midpoint  of  the  initial  pressure  rise.  Each  transducer 
record  is  also  used  to  estimate  the  time  for  the  "onset"  of  significant 
pressurization,  which  is  somewhat  subjective  but  not  unreasonable.  These  data 
appear  in  Fig.  14  along  with  the  location  of  the  projectile  face  after  impact; 
particle  velocity  is  nearly  constant  at  127  m/s  before  a  deceleration  begins  beyond 
200  ms  (presumeably  due  to  vigorous  combustion).  By  this  time,  however,  the 
leading  compaction  front  has  reached  the  downstream  end  of  the  tube.  Although 
the  tube  eventually  ruptured,  there  was  no  transition  to  detonation. 

Predictions  from  the  model  are  compared  to  the  PDC  data  in  Fig.  14. 
There  are  two  observations  here  which  seem  to  have  general  applicability.  The 
first  concerns  the  under-prediction  of  propagation  speed  of  the  leading  compaction 
front.  The  PDC  data  imply  a  wave  speed  of  494  m/s,  but  present  theory  predicts 
385  m/s  for  an  inert  wave  in  which  the  aggregate  maintains  the  equilibrium  stress 
state.  However,  if  this  compaction  wave  induces  reaction  accompanied  by  full  heat 
release  (1104  cal/g  for  TS-3659),  then  conversion  of  only  0.44%  of  the  solid  would 
create  enough  gas  pressure  to  stiffen  the  aggregate  and  allow  it  to  support  the 
wave  speed  of  494  m/s.  The  magnitude  of  gas  pressure  required  to  support  the 
increased  wave  speed  is  nearly  independent  of  an  important  ambiguity  concerning 
the  amount  of  energy  released  in  the  wave-induced  reaction.  Reaction-generated 
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Figure  14  -  Distance-time  data  for  PDC-80;  160  m/s  impact  on  60.296  TMD 
TS-3659  propellant.  Data  from  Ref.  15.  Tabular  values  are  model  predictions 
downstream  of  indicated  wave. 
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gas  pressure  is  essentially  a  function  of  the  product  of  the  energy  released  and 
the  amount  of  solid  converted  to  gaseous  products.  Little  evidence  is  available  to 
determine  the  amount  of  energy  release  accompanying  this  reaction.  However, 
liberation  of  the  full  heat  of  combustion  of  the  solid  phase  means  that  the  product 
gases  trapped  in  the  porosity  of  the  aggregate  will  be  near  the  equilibrium  flame 
temperature  (except  for  minute  amounts  of  reaction).  Unless  the  aggregate  has 
been  compressed  above  90%  TMD,  confined  gases  at  these  temperatures  should 
produce  visible  radiation  which  could  be  seen  by  a  high-speed  camera  looking 
through  a  Lexan  tube  confinement.  For  TS-3659  propellant,  Glancy  et  al.  [15] 
attempted  three  such  runs  at  impact  speeds  up  to  290  m/s  and  failed  to  see  any 
evidence  of  light  on  the  camera  records,  although  the  tube  ruptured.  This 
observation  implies  that  wave-induced  reaction  in  TS-3659  releases  only  a  fraction 
of  the  total  energy,  and  this  fraction  is  less  than  the  value  which  would  produce 
gas-phase  temperatures  of  say  1500°K.  As  an  example,  the  computational  results 
in  Fig.  14  show  that  a  compaction  wave  speed  of  494  m/s  could  be  supported  by 
reaction  of  1.32%  releasing  only  304  cal/g  (27.5%  of  the  total)  which  leads  to  a  gas 
temperature  of  1243°K. 

The  second  observation  concerns  the  implied  propagation  speed  and  strength 
of  the  second  compression  wave.  Differencing  the  two  data  points  determined  from 
the  transducer  records  as  the  "onset”  of  significant  pressurization  (see  Fig.  14) 
produces  a  wave  speed  of  624  m/s.  If,  in  the  model,  the  piston  velocity  is 
assumed  to  increase  by  1  m/s  at  95  ms,  this  anchors  a  very  weak  "shock"  wave 
which  propagates  into  the  aggregate  at  613  m/s  virtually  passing  through  the  data 
points.  The  implication  is  that  the  leading  edge  of  the  pressurization  wave 
propagates  near  the  local  sound  speed  in  the  compressed  material.  Since  the 
transducers  testify  to  a  rapidly  rising  pressure  field  downstream  of  wave  (2),  the 
quasi-steady  assumption  in  the  model  breaks  down  in  this  important  region.  Thus 
the  single  weak  wave  (2)  shown  in  Fig.  14  does  not  properly  represent  the 
influence  of  vigorous  combustion  which  begins  near  the  piston  face.  Successive 
pressure  increases  generated  by  transient  combustion  behind  wave  (2)  will  propagate 
at  increasing  speeds  toward  the  collision  point,  and  their  coalescence  may  influence 
the  "trigger"  mechanism. 

P DC-73  /  49.4%  TMD  WC-231  Propellant 


WC-231  is  a  double-base  ball  propellant  (nominally  25%  NG  in  NC)  whose 
grains  are  shaped  like  "fat  pancakes"  (spheres  which  have  been  rolled)  with  an 

average  diameter  of  790  Mm.  In  PDC-73  [15],  49.4%  TMD  WC-231  is  confined  in  a 
Lexan  tube  and  impacted  at  a  projectile  speed  of  220  m/s.  Transition  to 
detonation  did  not  occur  within  the  experimental  tube  length  (146  mm). 

A  simulation  is  shown  in  Fig.  15  along  with  the  PDC  data  [15].  If  the 

piston  velocity  is  equated  to  the  measured  particle  velocity  of  183  m/s,  the  model 
predicts  the  leading  compaction  front  (if  inert)  would  propagate  at  390  m/s;  the 

PDC  data  imply  approximately  430  m/s.  If  wave-induced  reaction  releases  the  full 
heat  of  combustion  (1150  cal/g,  assumed),  then  reaction  of  only  022%  of  the  solid 
would  support  the  observed  wave  speed.  The  luminosity  data  suggest  that  vigorous 
reaction  begins  near  the  piston  face  at  approximately  125  ms.  Assuming  a  1  m/s 
increase  in  piston  velocity  at  this  time  generates  a  weak  shock  wave  propagating 
at  736  m/s  which  forms  a  boundary  for  the  leading  edge  of  luminosity.  Although 
overly  simplistic,  if  collision  of  this  second  shock  wave  with  the  leading  compaction 
front  were  to  trigger  about  5%  reaction,  the  combined  shock  wave  would  propagate 
at  121  km/s  and  pass  through  the  experimental  data  for  strong  luminosity.  The 
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Figure  15  -  Distance-time  data  for  PDC-73;  220  m/s  impact  on  49.4%  TMD 
WC-231  propellant.  Data  from  Ref.  15.  Tabular  values  are  model  predictions 
downstream  of  indicated  wave. 
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predicted  magnitude  of  reaction  to  drive  the  combined  shock  wave  is  plausible,  but 
it  is  unlikely  that  collision  of  the  postulated  weak  second  shock  wave  with  the 
leading  wave  could  be  the  sole  trigger  for  the  increase  in  reactivity;  transient 
combustion  behind  the  second  wave  must  play  a  role.  The  fact  that  luminosity 
was  not  observed  directly  behind  the  leading  front  again  suggests  partial  energy 
release  in  the  wave-induced  reaction.  If  it  is  assumed  that  only  250  cal/g  are 
liberated  by  reaction  of  0.896  of  the  solid,  the  wave  pattern  in  Fig.  15  is  unchanged 
except  that  the  temperature  of  the  trapped  gases  reaches  only  1205°K.  In  both 
simulations  of  PDC-73,  the  leading  edge  of  luminosity  propagates  at  the  local 
sound  speed  in  the  aggregate  just  as  the  onset  of  pressurization  did  in  PDC-80. 

Partial  confirmation  of  the  reactive  wave  assumption  in  the  present  model 
may  be  provided  by  an  experimental  measurement  [15]  unique  to  PDC-73.  This 
run  included  two  infared  detectors  which  observed  the  aggregate  through  a  NaCl 
window  aperature  (50  mm  from  the  impacted  end)  in  the  wall  of  the  Lexan  tube. 
Although  the  IR  signals  were  not  calibrated  for  absolute  temperature,  their  ratio 
indicates  relative  temperature  [15].  First  response  from  the  two  detectors  occurred 
at  121  and  127  ns ,  while  the  maximum  ratio  (maximum  temperature)  occurred  at 
130  ns;  arrival  time  of  the  leading  compaction  wave  is  predicted  to  be  117  ns. 
This  strongly  suggests  that  the  compaction  wave  itself  initiated  reaction.  Just 
behind  the  wave  representing  the  leading  edge  of  luminosity,  the  IR  signals  indicate 
a  broad-based  participation  in  reaction  but  temperatures  remain  below  the  previous 
maximum  until  after  the  appearance  of  the  brightly  luminous  front  which  coincides 
with  the  combined  shock  wave  in  the  present  model. 

PDC-27  /  7396  TMD  Class  D  HMX 


Simulations  of  impact  on  Class  D  HMX  are  clouded  by  the  influence  of 

particle  size  distribution,  which  introduces  a  variability  in  quasi-static  compaction 
behavior  [38]  and  hence  wave  propagation  speed  in  the  present  study.  But  a 
simulation  of  one  of  Sandusky's  PDC  runs  will  be  used  here  to  illustrate  some 
consequences  of  the  multiple  wave  structure.  PDC-27  [14]  involved  Class  D  HMX 
loaded  at  7396  TMD  into  a  147  mm  long  Lexan  tube  and  impacted  at  267  m/s;  a 
transition  to  detonation  did  take  place  as  indicated  in  Fig.  16.  Since  luminosity  is 
coincident  with  the  leading  compaction  front,  the  simulation  of  wave-induced 

reaction  assumed  total  energy  release  and  found  that  reaction  of  approximately 

0.5%  is  required  to  support  the  observed  wave  speed  of  778  m/s.  If  collision 

triggers  reaction  of  3.2%,  the  combined  shock  wave  will  be  driven  at  2  km/s  and 
pass  through  the  next  four  data  points  for  strong  luminosity.  In  this  run,  the 
wave  continues  to  accelerate  and  quickly  achieves  the  detonation  velocity  of  7.46 
km/s. 


Up  to  this  point,  the  role  played  by  the  rearward  propagating  shock  wave 
(3)  and  the  contact  discontinuity  between  waves  (3)  &  (4)  has  been  rather  obscure. 
This  study  would  contend,  however,  that  the  density  differences  bounded  by  these 
waves  should  be  clearly  discernable  in  a  radiograph.  Fig.  17  is  a  plot  of  the 
spatial  distribution  of  mixture  density  (ratioed  to  initial  mixture  density)  at  35  #s, 
derived  from  the  computed  flow  field  for  PDC-27  shown  in  Fig.  16.  It  simulates 
what  a  flash  radiograph  would  record  if  triggered  (at  35  ns)  shortly  after  the 
collision.  There  are  several  interesting  features  here.  Density  of  the  reacting 
mixture  directly  behind  the  combined  shock  wave  (4)  traveling  at  2  km/s  is  not 
much  greater  than  the  quiescent  value.  Even  though  this  reacting  mixture  is 
supporting  a  stress  of  300  MPa,  the  solid  volume  fraction  is  virtually  unchanged 
because  the  high  pressure  gas-phase  combustion  products  can  carry  the  same  "load" 

(text  continues  on  page  38) 
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I^ure  16  -  Distance-time  data  for  PDC-27;  267  m/s  impact  on  7396  TMD  Class  D 
HMX.  Data  from  Ref.  14.  Tabular  values  are  model  predictions  downstream  of 
indicated  wave. 
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Figure  17  -  Simulated  radiograph  triggered  at  35  fis  in  PDC-27.  Numerical  values 
from  model  predictions  listed  in  Figure  16. 
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Figure  18  -  Simulated  radiograph  triggered  at  35  ns  in  ALTERED  PDC-27. 
Numerical  values  from  model  predictions  listed  in  Figure  16. 


as  the  solid  material.  The  high  density  region  ("plug")  is  material  trapped  behind 
the  slow  moving  contact  discontinuity.  The  highest  density  represent*  material 
which  has  been  compressed  by  shock  waves  (1),  (2)  &  (3).  The  "reverse  step"  in 
density  is  created  by  shock  wave  (3)  which  is  propagating  rearward  into  material 
already  subjected  to  (1)  &  (2).  This  density  pattern  is  accentuated  when  the 

collision  triggers  an  increased  level  of  reaction  in  the  combined  shock  wave  (4). 

For  demonstration  purposes,  the  flow  field  for  PDC-27  was  altered  by  increasing 
the  amount  of  reaction  assumed  to  be  triggered  by  the  collision.  Although  shock 

waves  (1)  &:  (2)  remain  unchanged,  the  case  labelled  ALTERED  PDC-27  (see  table 
under  Fig.  16)  predicts  a  combined  shock  wave  (4)  propagating  at  6  km/s.  The 

corresponding  simulated  radiograph  at  35  us  is  shown  in  Fig.  18.  Note  that 

mixture  density  behind  shock  wave  (4)  is  now  virtually  equal  to  the  quiescent 
value.  The  highest  density  in  the  "plug"  has  increased,  the  velocity  of  the  leading 
edge  (contact  discontinuity)  of  the  "plug"  has  slowed  to  less  than  100  m/s,  and  the 
reverse  step  is  more  pronounced  because  the  rearward  propagating  shock  wave  (3) 
is  stronger. 

PDC-27  did  not  happen  to  include  flash  radiography,  but  it  may  be 
appropriate  to  compare  the  above  prediction  to  a  result  obtained  by  McAfee  [10,11] 
in  64%  TMD  Class  A  HMX.  In  McAfee's  experiment,  combustion  of  ignitor 

material  at  one  end  in  a  closed  tube  creates  a  rising  pressure  field  against  a  steel 
cylinder  or  "piston"  which  is  then  driven  into  the  confined  Class  A  HMX. 

Confinement  is  a  fairly  thin- walled  maraging  steel  tube  which  allows  flash 
radiography  to  record  the  position  of  tracer  foils  seeded  into  the  original  bed; 

diagnostics  also  include  capped  and  ionization  pins.  The  X-T  plot  of  McAfee's  run 
B-9153  reproduced  here  as  Fig.  19a  shows  that  a  transition  to  detonation  (6.42 
km/s)  did  occur.  Radiograph  #2  was  taken  just  after  the  transition,  and  the 
results  are  reproduced  in  Fig.  19b.  Some  striking  similarities  are  evident  when 
comparing  this  experimental  result  with  the  simulated  radiograph  in  Fig.  18  for  the 
6  km/s  combined  shock  wave  (4).  If  it  can  be  assumed,  in  McAfee's  experiment, 
that  a  shock  wave  collision  triggered”  an  abrupt  increase  in  wave-induced  reaction, 
then  the  present  model  predictions  would  appear  to  be  compatible  with  his 
observations. 

DDT  SHOT  A266  /  59.8%  TMD  TS-3659 


At  the  present  time,  none  of  the  impact  experiments  have  produced  pressure 
transducer  records  which  offer  conclusive  evidence  for  the  presence  of  the  rearward 
propagating  wave.  Such  evidence  is  provided  by  a  DDT  experiment.  Bernecker  [8] 
has  reported  results  for  59.8%  TMD  TS-3659  confined  in  the  thick  wall  steel  tube 
which  was  instrumented  with  a  series  of  ionization  pins  as  well  as  several  wall- 
mounted  transducers  [8].  The  arrival  time  of  the  strong  reaction  front  can  be 
extracted  from  the  pin  data,  but  the  transducers  continue  to  record  a  time-history 
after  the  front  has  passed.  The  results  [8]  from  Shot  A266  are  reproduced  here  in 
Fig.  20,  in  terms  of  x*  and  t*  which  are  distance  and  time  referenced  to  the 
onset  of  detonation.  The  acceleration  of  the  strong  compaction  front  from  620 
m/s  to  2.15  km/s  is  accompanied  by  a  rearward  propagating  compression  wave 
(speed  of  ~1.7  km/s)  as  indicated  by  the  transducer  excursions  at  symbols  D,  E,  F 
&  G. 


The  present  theory  based  on  simple  motion  of  a  piston  driver  is  ill- 
equipped  to  simulate  a  DDT  experiment  which  begins  with  ignitor  combustion.  In 
particular,  no  experimental  value  of  particle  velocity  is  available  to  determine  the 
motion  of  a  hypothetical  "piston".  A  simulation  was  based  on  a  230  m/s  particle 
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velocity  and  0.2996  reaction  to  create  the  strong  compaction  front  at  630  m/s,  a  45 
Us  time  delay,  and  then  a  1  m/s  increase  in  particle  velocity  to  anchor  the  second 
wave.  Although  clearly  not  unique,  this  wave  structure  predicts  that  a  collision 
which  triggers  — 1096  reaction  would  produce  a  combined  shock  wave  (4)  traveling 
at  2.15  km/s.  Furthermore,  a  rearward  propagating  shock  wave  (3)  will  travel  at 
1.72  km/s  if  it  induces  an  additional  0.7296  reaction  in  the  already  burning  material 
compressed  by  waves  (1)  &  (2).  The  path  of  both  these  waves  lies  parallel  to  the 
experimental  data.  Comparison  of  the  transducer  time-history  to  the  predicted 
pressure  level  of  350  MPa  behind  shock  wave  (3)  indicates  that  theory  is  low  by 
approximately  a  factor  of  two,  which  may  not  be  surprising  when  considering  the 
assumptions  which  have  been  made.  Although  the  analysis  suggests  that  the 
rearward  wave  was  triggered  by  a  collision,  this  computational  example  does  not 
constitute  proof. 


Figure  20  -  Relative  distance-time  data  from  Bernecker's  [8]  DDT  Shot  A266  using 
59.8%  TMD  TS-3659  ball  propellant  (circles  =  probe  response,  squares  =  transducer 
excursion;  velocity  in  mm/#s). 


To  summarize  this  section,  various  laboratory  experiments  on  SDT  and  DDT 
in  granular  energetic  material  have  observed  abrupt  increases  in  the  speed  of 
propagation  of  compaction  waves  and  the  leading  edge  of  luminosity.  The  DDT 
experiment,  in  at  least  one  case,  has  evidence  of  a  strong  rearward  propagating 
wave.  A  possible  explanation  is  provided  by  a  simple  quasi-steady  model  which 
envisions  an  increase  in  reactivity  triggered  by  the  collision  of  two  reactive 
compaction/shock  waves.  Simulation  of  experimental  data  shows  that 

(a)  The  leading  compaction  wave  created  by  modest  impact  on  confined 
granular  energetic  material  is  not  necessarily  inert  The  experimental  propagation 
speeds  which  exceed  the  predicted  inert  wave  speed  could  be  supported  by  wave- 
induced  reaction  of  a  few  tenths  of  one  percent  of  the  solid  material,  if  the 
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reaction  releases  the  total  heat  of  combustion.  If  wave-induced  reaction  is 
accompanied  by  partial  heat  release  (highly  probable  for  single  and  double  base 
propellants),  then  the  magnitude  of  reactivity  is  near  one  percent. 

(b)  After  a  time  delay,  the  onset  of  vigorous  combustion  of  material  adjacent 
to  the  piston  face  will  drive  a  second  compressive  wave  system  into  the  aggregate 
which  has  been  formed  by  the  initial  compaction  wave.  The  leading  edge  of  this 
second  wave  system  propagates  at  a  speed  near  the  local  sound  velocity. 

(c)  The  analysis  predicts  that  the  abrupt  increases  in  propagation  speed 
observed  experimentally  would  require  the  collision  to  induce  reaction  in  the  range 
of  2%  -  10%.  Transient  combustion  (neglected  in  this  model)  behind  the  leading 
edge  of  the  secondary  wave  system  must  play  a  role  in  the  "trigger”  mechanism. 

(d)  The  wave  pattern  produced  by  the  collision  includes  both  a  combined 
compaction/shock  wave  exhibiting  an  abrupt  increase  in  speed,  and  a  rearward 
wave  which  propagates  back  into  material  compressed  by  the  first  two  waves.  A 
simulated  radiograph  taken  shortly  after  the  collision  predicts: 

(1)  The  density  behind  the  combined  shock  wave  may  appear  nearly  equal 
to  the  quiescent  value.  Wave-induced  reaction  generates  enough  gas  pressure  to 
counteract  (resist)  the  normal  occurrence  of  bed  compaction  caused  by  a  large 
value  of  mixture  stress. 

(2)  A  high  density  ("plug”)  region  is  trapped  behind  the  slow-moving  contact 
discontinuity. 

(3)  The  spatial  distribution  of  the  high  density  region  will  exhibit  a  "reverse 
step"  as  a  direct  result  of  the  rearward  propagating  shock  wave. 

It  appears  from  these  simulations  that  the  appearance  of  a  rearward 
propagating  wave  may  be  a  direct  indicator  that  a  "trigger"  mechanism  is 
responsible  for  the  abrupt  increase  in  wave  speed.  The  lack  of  any  rearward 
propagating  wave  may  signify  that  the  transition  mechanism  is  controlled  by  nearly 
complete  combustion  which  "drives"  from  behind  until  it  catches  the  leading 
compaction  front,  which  is  essentially  the  mechanism  deduced  fifteen  years  ago  by 
Bernecker  &  Price  [4]. 
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VI.  TRANSIENT  REACTIVE  SHOCK  WAVE  MODEL  [36] 

A  basic  objective  here  is  to  create  a  more  realistic  description  of  the 
transient  events  seen  in  the  PDC  experiment  by  removing  many  of  the  restrictive 
assumptions  adopted  in  the  previous  three  models.  This  one-dimensional  model  is 
designed  to  simulate  the  configuration  shown  in  Fig.  1.  As  before,  the  mixture  of 
compressible  granular  solids  and  gas  deforms  in  uniaxial  strain.  It  is  also  assumed 
that  the  mixture  responds  to  a  change  in  density  by  instantaneously  adjusting  to 
the  equilibrium  stress  state  defined  in  Section  II,  Le.,  rate-dependent  resistance  to 
porosity  change  is  neglected.  However,  this  model  does  account  for 


(a)  deviator  stress  components  in  the  formulation  of  the  mixture  stress  tensor, 

(b)  a  traction  force  (friction)  between  the  sliding  aggregate  and  the  confining 
boundary, 

(c)  a  drag  force  between  the  gas  and  solid  particles  (gas  and  solid  phases  have 
different  velocities), 

(d)  heat  transfer  between  the  two  phases, 

(e)  ignition  and  transient  combustion  of  the  solid  phase,  and 

(f)  reaction  of  a  gas-phase  intermediate  species  to  a  final  combustion  product 

Although  the  impact-generated  compaction  wave  is  modeled  as  a  shock  wave,  no 
propagation  speed  is  constrained  to  be  quasi-steady.  Combustion-generated 
compression  wave  systems  can  appear  at  any  time  in  the  transient  flow  field. 

The  conversion  of  energetic  solid  material  to  final  gaseous  combustion 
products  is  a  complex  process  which  is  not  well  understood.  Although  solid 

propellant  combustion  can  often  be  represented  with  a  pressure-dependent 
prescription  for  surface  regression  rate,  this  implicity  assumes  that  the  gas-phase 
flame  zone  is  quasi-steady  and  remains  anchored  to  the  interface.  Considerable 

additional  detail  is  required  for  a  description  of  ignition  and  the  transient  events 
leading  up  to  an  "established"  flame  zone  near  the  solid  surface.  Even  if  the 

complete  decomposition  scheme  and  associated  multi-species  reaction  parameters 
were  known,  this  level  of  detail  is  far  beyond  the  scope  of  the  present  effort  A 
first  step,  however,  is  the  simple  two  species  reaction  scheme  sketched  in  Fig.  21. 
It  is  assumed  that  "Rj"  is  a  reactive  intermediate  gas-phase  species  which  is  the 
result  of  a  partial  or L  incomplete  decompostion  of  the  energetic  solid.  Then,  "R," 
will  liberate  the  remaining  energy  in  a  gas-phase  reaction  to  form  "R,",  the  final 
flame  product  The  present  model  envisions  a  staged  combustion  process  which 
changes  at  a  "transition"  value  of  gas  pressure  defined  as  P.  (assumed  constant 

here).  When  the  local  gas  pressure  is  less  than  ,  solid-pnase  decomposition  is 
incomplete  and  releases  the  reactive  intermediate  species  R.  into  the  surrounding 
gas  mixture.  When  the  local  gas  pressure  exceeds  P,,  solid-phase  decomposition 
goes  to  completion  and  releases  the  final  product  gas  R«.  The  gas-phase  reaction 
which  converts  R.  to  Rg  can  occur  at  any  time  or  location.  The  precedent  for 
such  a  scheme  can  be c  found  in  the  two-stage  flame  zone  of  a  double-base 
propellant,  and  in  the  ignition  and  flame  spreading  problems  [51]  observed  with 
certain  LOVA  nitramine  propellants  used  in  guns.  Price  et.  al.  [23,24]  have  been 
pursuing  similar  ideas  for  several  years. 
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REACTION  /  COMBUSTION  SCHEME 

R[  =  Reactive  Intermediate  at  Energy  Level  Aegr 
R£  =  Final  Combustion  Product  at  Energy  Level  eg 
In  Gas  Phase:  Rj  "" *■  R2 

For  Gas  Pressure  <  Transition  Pressure: 

n  Gas-Phase 
R1  - -  R2 


For  Gas  Pressure  >  Transition  Pressure: 


Figure  21  -  Schematic  of  combustion  process  in  modeL 


The  model  reaction  scheme  is  an  extreme  simplification,  but  it  may  be 
useful  in  demonstrating  potential  interactions.  The  impact-generated  compaction 
wave  begins  the  process  by  inducing  reaction  of  a  small  amount  of  energetic  solid 
within  the  wave  front  itself  (see  discussion  in  Section  IV).  It  is  assumed  that  this 
wave-induced  reaction  produces  the  reactive  intermediate  gas  species,  R^,  which  is 
trapped  in  the  porosity  of  the  compressed  aggregate.  This  trapped  gas1  must  wait 
through  an  induction  time  or  delay  time,  t.  .,  before  the  reaction  which  converts 
R.  to  R«  can  begin;  the  idea  of  a  delay  time  was  taken  from  the  work  of  Baer 
and  Nunziato  [20].  Once  the  gas-phase  reaction  is  underway,  the  pressure  and 
temperature  of  the  gas-phase  mixture  will  increase.  This  increases  the  heat 
transfer  rate  to  the  solid  phase,  which  raises  the  solid-phase  surface  temperature 
(predicted  with  the  scheme  devised  by  Baer  [48]).  Solid-phase  combustion  begins 
when  the  computed  surface  temperature  exceeds  an  ignition  value  (typically,  473°K 
here).  Combustion  of  the  solid  phase  proceeds  at  a  prescribed  pressure-dependent 
rate  producing  intermediate  species  Rj  when  P  <  P^.  and  then  R2  above  the 
transition  pressure.  Given  an  initial  value  of  uelay  tune,  t^0 ,  the  local  value 
follows  from  a  solution  of, 


D(tde,)/Dtg  =  -<Pg/Pg0>2- 


(15) 


The  global  gas-phase  reaction,  R^  — > 


proceeds  at  the  finite  rate,  m 


gr’ 


where 


for  t,  ,  >  0.0  m _ =  0.0 

del  gr 

for  t ,  ,  <  0.0  m_  =  A  Y.  p  exp(  -E  /  T  ) 
del  gr  gr  I  rg  r'  gr  g 


(16) 


Since  this  is  not  an  elementary  reaction,  the  rate  parameters  A  and  E  must  be 
estimated.  Keeping  in  mind  that  the  mass  generation  rate,  m,  from  ^olid-phase 
combustion/pyrolysis  produces  R,  when  P  <  P^r  and  R£  otherwise,  the  time 
rate-of-change  of  Yj  can  be  written,  ° 
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D(Yj)/Dtg  =  [  m(when  Pg<Ptr)  -  mgr  -  Yj  m  J  for  <  0.0  (17) 

where  the  first  term  on  the  RHS  vanishes  when  the  solid  combustion  process  goes 
directly  to  completion  (P  >  Ptr). 

The  balance  equations  are  written  in  a  translating  coordinate  system  where 

V  =  V  (t)  =  local  velocity  of  axial  coordinate, 
c  c 

Relative  velocities  for  each  phase  are  then 

Ug  =  ug  -  Vr  =  solid-phase  velocity  relative  to  local  coordinate  system, 


and 


Ug  =  ug  -  Vc  =  gas-phase  velocity  relative  to  local  coordinate  system. 


Balance  of  Mass  /  Solid  Phase: 

DPS 


D  t 


+  Ps  US-  +  PS 


D  t 


O 

-  m 


(18) 


Balance  of  Mass  /  Gas  Phase: 

DP, 


M  Dt 


+  -°g  ugz  -  P& 


Df. 

Dt 


=  +  m 


(19) 


Balance  of  Momentum  /  Solid  Phase: 

Du. 


D  t 


+  °SZ  +  £SZ  =  fd  -  f. 


(20) 


Balance  of  Momentum  /  Gas  Phase: 


Du, 

'*  V  Dt. 


+  £bP.  = 

8  8z 


-  fd  -  m  (ug  -  us) 


(21) 


Balance  of  Energy  /  Solid  Phase: 
D  e. 


D  t0 


+  esffs\+ 


D ^ 

D  t„ 


O 

m 


=  q -(<%-£»)- 

"s 


(22) 
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Balance  of  Energy  /  Gas  Phase: 


'«  V  D  t 


D  £, 


+  £„  P„  U„  —  P 


«  «  *Z  *  \  Dt 


+  fd(ug-us)  -q  +  m  (es-eg)  + 


[ 


(CTs  -  fin)  (us  - 

- - -  + 


\2  “I 

US  -  U«) 

2 


where  —  -  -  =  (  ).  +  U  (  ), 

Dt  Wt  sV  h 


TT"()t  +  u«  (>■ 


(23) 


The  above  equation  system  is  incomplete  without  a  relationship  between  porosity 
and  the  state  of  stress  in  the  mixture;  in  the  present  analysis,  this  rate- 
independent  relationship  is  denoted  the  "equilibrium  stress  state"  and  is  described  in 
Section  II  [see  Eq.  8). 

A  deliberate  attempt  was  made  to  preserve  the  integrity  of  the  equation 
system  by  deriving  a  numerical  solution  procedure  based  on  the  method-of- 
characteristics  (MOC).  A  MOC  solution  holds  the  promise  of  faithfully  reproducing 
the  local  wave  motion  and  propagating  shock  waves  as  discontinuities  without 
spurious  numerical  distortion.  The  disadvantages  include  a  complicated  derivation 
and  a  non-trivial  coding  exercise  to  obtain  the  solution.  These  disadvantages  are 
accentuated  in  two-phase  flow.  Single-phase  flow  yields  three  characteristic 
directions  and  their  associated  compatibility  conditions,  all  of  which  can  be  solved 
explicitly  although  the  procedure  involves  iteration.  The  current  system  of  two- 
phase  flow  equations  which  are  constrained  to  the  equilibrium  stress  state  is  more 
difficult  to  solve.  The  system  consists  of  a  mass,  momentum,  and  energy  equation 
for  each  phase  (Eqs.  18-23),  along  with  a  statement  of  the  equilibrium  stress  state 
(Eq.  8). 


Given: 


P*  (*.)  =  <%  -  p,  ;  <7,  =  p»  K' e»)  /  ^  (£») 


d 


de. 
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Define: 


a«2  =  (pe)„/(p</0(pl)e< 

a..2s[(p*)„+(as/ps2)(p.)..]  /* 

A,.2-al .2  -  (^  /  W)  (P>)e 

3 

Pga«2  psais2  PsAis2  °sr 

8  £*  £s  A  e3  f 

C2S  («A  +  P*  )  /  {*g  +  «K  +  Ps) 

Cy  =  [°g  +  ^8  +  (^s  /  £s)]  /  (ag  +  aA  +  ^8  ) 

C?  S  (^s  +  /  £s)  /  (a£  +  aA  +  ) 

The  six  distinct  characteristic  directions  follow  from  the  roots  of 

ws  W  (  ws2  W  2  -  C~  a  2  wa2  -  C~  a,  2  W  2  +  C~  a  2  als2  )  =  0 

sg\sg  xgs  ylsg  zglsy 

dz  dz 

where  W„  = - U  ;  W„  = - U 

s  dt  8  *  dt  g 


(24) 


The  two  trivial  roots,  W  =  0  and  W  =  0,  represent  the  respective  streamlines. 
Unfortunately,  the  fourth-order  expression  in  brackets  does  not  factor  explicitly, 
and  the  four  "acoustic"  directions, 


W  +  ,  W" 
g  g 


are  the  four  real  roots  which  must  be  found  numerically.  The  relationship  of 

these  roots  is  sketched  in  Fig.  22  below,  where  7-1  represents  the  solid-phase 
streamline,  6-1  is  the  gas-phase  streamline,  2-1  and  5-1  are  the  right-running  and 

left-running  solid-phase  acoustic  directions,  respectively,  and  3-1  and  4-1  are  the 

right-running  and  left-running  gas-phase  acoustic  directions,  respectively.  The 

compatibility  conditions  face  a  similar  complication.  Both  energy  equations  can  be 
written  in  characteristic  form  along  their  respective  streamlines,  i.e^ 


P8  =  P8(°8)  and  pg  =  pg(Pg). 

However,  the  remaining  four  dependent  variables  (ugl  u  ,  og  and  eg,  here)  are 
coupled  in  the  matrix  equation,  8 
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M-O-C  Solution  Procedure 


Figure  22  -  Schematic  representation  of  method-of-characteristics  solution  procedure. 


Fs(V-«t2)-f'«K2-als2)  =  0 
where 

Fs-P,als2(u3-HU3)  +Ws[(^-H„s)  +cA(?5-Hes)] 

Fg  E  P.  ag2  (“g  ~  H»e)  +  Wg  [  (*.  -  H„.)  -  («.  +  A)  (i.  -  Hts)  ] 

and  C)-ft() 


(25) 


With  reference  to  Fig.  22,  the  solution  is  found  in  the  following  seqence.  Using  an 
estimate  for  all  six  characteristic  directions,  Eq.  25  is  inverted  to  give  values  for 
ug,  u  ,  a  and  €g  at  the  new  point  "1”.  The  stress  state  equation  immediately 
specifies  F  .  Thestwo  energy  equations  yield  the  density  values,  and  the  equations 
of  state  provide  the  internal  energies.  Then,  the  characteristic  directions  must  be 
re-solved  (Eq.  24)  and  the  whole  process  repeated  until  convergence  is  obtained. 
Although  this  procedure  is  somewhat  involved,  flow  fields  for  the  problems 
discussed  below  require  only  1  to  2  minutes  on  the  BRL  Cray-2. 

Within  the  constraints  of  the  present  review,  it  is  not  possible  to  include  a 
discussion  of  all  the  computations  made  with  this  model  Instead,  this  section  will 
focus  on  three  important  areas  in  the  simulation  of  NSWCs  PDC  experiment:  (A) 
the  influence  of  +he  projectile  path,  (B)  simulation  accuracy  when  the  granular 
material  is  inert,  and  (C)  the  influence  of  the  energy  release  mechanism.  In  (B), 
the  examples  are  based  on  PDC-M34  and  PDC-74  which  involve  the  inert  material 
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melamine  (46-56  Mm  diameter).  In  (A)  and  (C),  the  examples  are  based  on  the 
experimental  runs  PDC-80  and  PDC-81,  both  of  which  involve  Lexan  projectile 
impact  on  approximately  6096  TMD  TS-3659  ball  propellant  (spherical  grains,  434 
Mm  diameter)  which  has  a  nominal  composition  of  2196  NG  in  NC  (estimated  heat 
of  explosion  is  1104  cal/g).  Additional  details  about  the  experiments  can  be  found 
in  Refs.  15  and  42;  both  PDC-80  and  81  were  part  of  the  JANNAF  Modeling 
Workshop  summarized  in  Ref.  32. 

A.  Influence  of  Projectile  Path  /  PDC  Experiment 

Typical  low-range  response  from  the  two  wall-mounted  pressure  transducers 
can  be  found  in  PDC-81  (see  Fig.  23b;  gage  #1  is  located  at  38.1  mm,  gage  #2  is 
at  762  mm,  and  tube  length  is  146.8  mm).  After  the  impact-generated  compaction 
front  passes  each  gage  location,  the  response  shows  a  nearly  linear  "decay"  from 
the  initial  maximum  value  until  the  beginning  of  rapid  pressurization  presumably 
due  to  significant  reaction.  The  strength  of  this  linear  decay  seems  to  increase 
with  projectile  speed  at  impact  None  of  the  models  at  the  JANNAF  Modeling 
Workshop  [32]  predicted  this  decay,  and  the  behavior  went  unexplained.  Predictions 
of  a  runaway  event  depend  critically  on,  and  are  extemely  sensitive  to,  the 
competition  between  sources  of  heat  generation  (reactions)  and  heat  loss  mechanisms 
which  occurs  in  this  region.  A  successful  model  must  be  able  to  predict  this 
stress  state. 

The  results  displayed  in  Figs.  23-25  are  simulations  of  PDC-81  compared  to 
experimental  data  for  compaction  wave  locus,  projectile  front  locus,  and  the  wall- 
mounted  gage  records.  The  simulations  assume  no  gas-phase  reaction,  and  no 
solid-phase  combustion  other  than  a  small  amount  of  reaction  induced  by  the 
compaction  wave  itself.  Note  that  the  experimental  data  for  projectile  face 
location  (squares,  in  Fig.  23a)  form  the  basis  for  the  time-dependent  boundary 
condition  at  the  impacted  boundary.  Assuming  the  projectile  boundary  to  have  a 
constant  velocity  of  192  m/s  produces  the  results  shown  in  Fig.  23.  Recall  that 
the  present  theory  assumes  an  instantaneous  adjustment  to  the  equilibrium  stress 
state  and  predicts  compaction  waves  to  be  infinitesimally-thin  shock  waves.  Since 
true  wave  thickness  is  not  simulated,  the  best  possible  comparison  is  for  the 
predicted  pressure  jump  behind  the  shock  wave  to  occur  at  the  mid-point  of  the 
initial  rise  reported  by  the  experimental  gage.  As  seen  in  Fig.  23,  the  locus  of  the 
predicted  compaction  wave  is  close  to  the  NSWC  microwave  interferometry  data, 
and  the  timing  of  the  initial  response  from  the  simulated  wall  gages  is  close  to 
the  mid-point  of  the  actual  gage  rise  time.  Similar  to  many  other  cases,  the 

predicted  magnitude  of  the  initial  maximum  value  is  also  close  to  the  experimental 
value,  implying  that  these  gages  are  responding  to  the  radial  component  of  the 
equilibrium  mixture  stress  state.  The  predicted  gage  response  shows  no  evidence  of 
a  "decay”;  instead  there  is  a  slight  increase  as  the  result  of  wall  friction  operating 
on  a  continuously  increasing  amount  of  aggregate  which  is  sliding  past  the  lateral 
confining  boundary.  Of  course,  with  the  assumption  of  a  constant  projectile 

velocity,  this  computation  did  not  accurately  simulate  the  projectile  face  location 
beyond,  say,  50  ms. 

If  it  is  assumed  that  the  projectile  velocity,  near  40  ms,  slows  down 

smoothly  from  192  m/s  to  173  m/s,  the  predicted  projectile  locus  ("dual-linear") 
closely  matches  the  experimental  data  as  illustrated  in  Fig.  24a.  Now,  however,  the 
computed  compaction  wave  locus  falls  increasingly  behind  the  microwave  data,  and 
the  computed  response  at  gage  #2  is  late  compared  to  the  experimental  gage 

record  (Fig.  24b).  The  expansion  wave  created  by  the  slowing  projectile  has 

(text  continues  on  page  54) 
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Figure  23  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-81: 
237  m/s  impact  of  Lexan  projectile  on  60.196  TMD  TS-3659  ball  propellant.  Model 
computation  assumes  projectile  interface  speed  =  192  m/s  (constant),  and  no 
gas-phase  reaction  or  propellant  ignition. 


(SdW)  ssa-qs  rBTP^H  ampcijn 


Figure  23  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-81: 
237  m/s  impact  of  Lexan  projectile  on  60.1%  TMD  TS-3659  ball  propellant.  Model 
computation  assumes  projectile  interface  speed  =  192  m/s  (constant),  and  no 
gas-phase  reaction  or  propellant  ignition. 
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Figure  25  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-81: 
237  m/s  impact  of  Lexan  projectile  on  60.196  TMD  TS-3659  ball  propellant.  Model 
computation  based  on  "quadratic"  projectile  path;  projectile  interface  speed  is  a 
quadratic  function  of  time.  Computation  assumes  no  gas-phase  reaction  or 
propellant  ignition. 
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Figure  25  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-81: 
237  m/s  impact  of  Lexan  projectile  on  60.1%  TMD  TS-3659  ball  propellant  Model 
computation  based  on  "quadratic"  projectile  path;  projectile  interface  speed  is  a 
quadratic  function  of  time.  Computation  assumes  no  gas-phase  reaction  or 
propellant  ignition. 


lowered  the  magnitude  of  both  gage  responses  (Fig.  24b),  but  at  best  the  computed 
response  represents  an  "average"  of  the  experimental  data.  The  linear  decay  region 
is  still  not  predicted. 

Out  of  frustration,  this  investigation  began  to  examine  the  problem  from 
another  point  of  view.  With  Ref.  26b  and  the  discussions  at  the  JANNAF 
Modeling  Workshop  [32]  as  background,  it  was  postulated  that  substantial 
deformation  of  the  front  face  of  the  Lexan  projectile  in  the  PDC  experiment  may 
indeed  be  important.  This  deformation  might  not  be  reflected  in  the  position  of 
scribe  lines  on  the  body  of  the  projectile,  which  is  the  basis  of  the  experimental 

data  for  projectile  locus.  Assume  for  the  moment  that  the  microwave  data  for 
compaction  wave  locus  is  more  accurate  than  the  data  given  for  projectile  locus. 
On  close  examination,  the  microwave  data  for  compaction  wave  locus  actually  show 
a  curved  path.  If  it  can  be  assumed  that  the  compaction  (shock)  wave  locus  up 
to  a  time  of,  say,  120  us  is  not  substantially  influenced  by  reaction,  then  the 

effective  projectile  path  must  also  have  been  curved  (not  two  straight  lines).  This 
study  posed  the  inverse  problem:  attempt  to  determine  the  projectile  path  which 
will  create  the  experimental  compaction  wave  locus.  One  possibility  (not  a  unique 
answer)  is  shown  by  the  results  in  Fig.  25.  As  can  be  seen,  a  slight  curvature  to 
the  projectile  path  (projectile  velocity  is  quadratic  in  time)  leads  to  a  computed 
compaction  wave  locus  which  is  a  close  match  with  the  NSWC  microwave  data 
(Fig.  25a).  Furthermore,  the  predicted  gage  response  DOES  reproduce  the  linear 
decay  behavior  seen  by  the  wall-mounted  pressure  transducers  (Fig.  25b).  On  the 

basis  of  these  examples,  this  investigation  feels  with  reasonable  certainty  that  the 

front  face  of  the  Lexan  projectile  in  the  PDC  experiments  was  deforming  under 
impact,  and  that  the  true  locus  was  a  curved  path.  The  important  implication  is 
that  the  true  path  is  not  accurately  known.  As  discussed  below,  this  has 
important  consequences  in  attempting  to  verify  reaction  mechanisms  when  the  basis 
of  comparison  is  the  wall-mounted  gage  records. 

B.  Simulation  of  PDC-M34  /  65%  TMD  Melamine  (Inert  Material) 


Before  addressing  the  complexities  created  by  impact-induced  reaction,  it 
would  be  interesting  to  evaluate  how  well  the  theory  can  simulate  the  NSWC  PDC 
experiment  when  the  confined  granular  material  is  inert.  A  recent  run  denoted 
PDC-M34  [53]  (see  Fig.  26  below)  involved  65%  TMD  melamine  (Eastman  1540) 
which  has  a  particle  size  distribution  between  46  -  56  Mm.  In  addition  to  the 
fact  that  melamine  should  behave  as  an  inert,  PDC-M34  employed  an  aluminum 

projectile  (206  m/s  initial  speed)  which  should  not  deform  under  impact.  This 

removes  the  uncertainties  in  projectile  locus  caused  by  deformation  of  a  Lexan 

projectile  as  discussed  in  Section  VI-A.  Furthermore,  confinement  in  PDC-M34 

was  provided  by  a  147  mm  long  thick-walled  steel  tube  which  effectively  eliminates 
radial  expansion  of  the  tube  wall.  As  in  the  earlier  PDC  experiments,  the  locus 
of  the  projectile  face  was  determined  from  the  motion  of  scribe  lines  on  the 

projectile  body  as  recorded  by  a  high-speed  camera  looking  through  a  slit  in  the 

tube  wall.  Microwave  interferometry  through  the  downstream  end  of  the  tube  was 
used  to  determine  the  location  of  the  strong  compaction  front.  The  time-history 
of  the  stress  field  was  monitored  at  three  fixed  locations  by  the  wall-mounted 
pressure  transducers  noted  in  Fig.  26. 

Before  the  transient  model  can  predict  wave  behavior,  it  is  necessary  to 

evaluate  the  equilibrium  stress  state  of  the  granular  aggregate.  This  calibration  is 

done  against  data  obtained  in  the  quasi-static  compaction  experiment.  Nearly  ten 
years  ago,  Elban  et.  al.  [52]  determined  the  three  data  sets  plotted  in  Fig.  27a  as 
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NSWC  Shot  PDC-M34  /  65%  TMD  Melamine 


Wall  Pressure  (Stress)  Time  History 
127.0  mm - 


-76.4  mm 


;38.2  mm  : 


ME 


Compacted 

Aggregate 


65%  TMD 
Melamine 


■ 

/ 


Compaction  Wave 


DIAGNOSTICS: 


Microwave 

Interferometry 


(a)  Microwave  Interferometry 

(b)  3  Wall-Mounted  Pressure  Transducers 

Figure  26  -  Schematic  of  NSWC  piston-driven-compaction  shot  PDC-M34  [53]. 


representative  of  65%  TMD  melamine.  The  solid  line  shown  in  Fig.  27a  is  a  "best 
fit"  approximation  generated  by  the  present  theory  when  r,=1.0  (kpsi),  Pj=0.0, 
B2=4.0  and  P2=0.7  (same  as  in  Fig.  A-l,  p.  67  in  Ref.  34).  in  the  discussion  to 
follow,  this  representation  will  be  referred  to  as  "regular"  melamine.  The 
equilibrium  shock  wave  theory  developed  in  Ref.  34  uses  the  equilibrium  stress 
state  to  predict  values  of  compaction  wave  speed  as  a  function  of  particle  velocity. 
These  predictions  for  regular  melamine,  shown  as  the  dashed  line  in  Fig.  28,  are 
compared  with  values  of  compaction  wave  speed  from  two  earlier  PDC  experiments 
(13]  (PDC-41A  and  41B)  in  65%  TMD  melamine.  The  comparison  is  well  within  the 
range  attributable  to  the  particle  size  distribution  in  melamine  as  illustrated  by  the 
quasi-static  compaction  data  sets  in  Fig.  27a.  Thus,  the  analysis  of  Ref.  34  which 
assumes  that  the  granular  aggregate  instantaneously  adjusts  to  the  equilibrium 
stress  state  appears  to  predict  the  correct  behavior  of  compaction  wave  speed  in 
melamine. 

The  PDC  experiment  of  interest  here  is  PDC-M34  [53]  which  also  involved 
65%  TMD  melamine  but  was  run  very  recently  (Dec  1989).  The  value  of 
compaction  wave  speed  (from  the  microwave  data)  just  after  impact  is  plotted  in 
Fig.  28  which  shows  that  this  value  is  located  considerably  below  the  prediction  for 
an  inert  compaction  wave.  Now  if  dynamic  resistance  to  porosity  change  at  this 
particle  velocity  has  become  important,  the  observed  value  of  wave  speed  should  lie 
above  (not  below)  the  theoretical  prediction  which  has  ignored  it  An  explanation 
may  be  related  to  the  age  of  the  material;  ten-year  old  melamine  may  have 
absorbed  moisture  and  aged  such  that  the  granular  aggregate  now  offers  less 
resistance  (is  "softer")  to  compaction  forces  than  it  did  in  Elban's  [52]  original 
measurements.  Exploring  this  hypothesis,  the  present  study  arbitrarily  adjusted  one 
of  the  parameters  (Bn=1.75)  controlling  the  strength  of  the  equilibrium  stress  state 
to  obtain  the  solid  Tine  shown  in  Fig.  27b  denoted  as  "soft"  melamine.  Again 
using  the  equilibrium  shock  wave  theory  [34],  predictions  for  compaction  wave  speed 
in  "soft"  melamine  are  shown  as  the  solid  line  in  Fig.  28;  the  intent  was  to  get 

(text  continues  on  page  62) 
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*  Quasi— Static  Compaction  *  65%  TUD  (r«f)  Halmwlna  *  Data:  Bban  •  N9VC 
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*  85%  TMD  Melamine  (soft)  *  INERT  SHOCK  WAVE  *  Data:  Sandusky  •  NSWC  * 
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interface  speed  (particle  velocity)  for  melamine  initially  at  6596  TMD.  PDC  data 
(41 A  and  41B)  from  Sandusky  and  Liddiard  [13];  dashed  line  is  prediction  based  on 
"regular”  melamine.  PDC  data  (M34)  from  Sandusky  and  Glancy  [53];  solid  line  is 
prediction  based  on  "soft”  melamine. 


*  INERT  SHOCK  WAVE  *  85%  TMD  (soft)  Melamine  *  Data:  Sandusky  •  NSWC 
1600.0  i  i  i  i  l  — '  i  . 1  I  l  I 


(s/m)  nsKHS  3AYM.  NOLIOVdHOD 


59 


Figure  29  -  Using  analysis  of  Ref.  34,  predicted  compaction  wave  speed  vs. 
interface  speed  (particle  velocity)  for  melamine  initially  at  8596  TMD.  PDC  data 
(74)  from  Sandusky  and  Glancy  [53];  solid  line  is  prediction  based  on  "soft" 
melamine. 
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Figure  30  -  Comparison  of  model  predictions  with  NSWC  data  [53]  from  PDC-M34: 
206  m/s  impact  of  aluminum  projectile  on  6596  TMD  melamine.  Simulation  ignores 
wall  friction  and  assumes  melamine  is  inert 
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Figure  30  -  Comparison  of  model  predictions  with  NSWC  data  [53]  from  PDC-M34: 
206  m/s  impact  of  aluminum  projectile  on  65%  TMD  melamine.  Simulation  ignores 
wall  friction  and  assumes  melamine  is  inert. 


close  to  the  experimental  value  for  PDC-M34,  not  to  match  it  exactly.  Another 
recent  PDC  run  in  melamine  may  support  the  "softening"  hypothesis.  During 
check-out  of  a  remotely  operated  quasi-static  compaction  apparatus  at  NSWC, 
melamine  was  compacted  to  85%  TMD.  The  PDC  run  denoted  PDC-74  [53] 
involved  a  217  m/s  impact  on  this  85%  TMD  aggregate.  Since  there  is  no 
quasi-static  compaction  data  available  for  melamine  initially  at  85%  TMD,  the 
present  theory  was  used  to  estimate  the  behavior,  and  then  this  equilibrium  stress 
state  was  supplied  to  the  equilibrium  shock  wave  theory  [34].  The  resultant 
predictions  are  shown  in  Fig.  29  along  with  the  experimental  value  of  compaction 
wave  speed  from  PDC-74.  The  comparison  is  quite  good,  suggesting  that  the 
assumed  "soft"  melamine  properties  do  describe  the  behavior  of  ten-year  old 
melamine  used  recently  in  the  PDC  experiment. 

The  simulation  of  PDC-M34  is  based  on  the  assumption  that  the  equilibrium 
stress  state  is  described  by  "soft"  melamine.  The  time-dependent  projectile-face 
boundary  condition  follows  from  a  quadratic  function  of  time  which  was  fit  to  the 
NSWC  projectile  path  data  (see  Fig.  30a).  The  predicted  compaction  wave  locus 
(solid  line  in  Fig.  30a)  is  slightly  ahead  of  the  NSWC  microwave  data,  but  this 

calculated  compaction  wave  locus  appears  to  mimic  the  observed  curvature.  The 
predicted  transducer  response  (radial  component  cf  mixture  stress)  is  compared  to 
the  experimental  data  for  all  three  gage  locations  in  Fig.  30b.  The  computed 
values  of  initial  response  seem  to  be  close  to  the  transducer  data  in  all  three 

cases.  Note  that  the  quadratic  path  of  the  aluminum  projectile  has  produced  a 
linear  "decay"  behavior  in  each  computed  gage  response  which  is  quite  similar  to 
that  discussed  in  Section  VI-A  above.  However,  there  are  significant  deviations 
after  the  initial  response.  Particularly  at  gage  #1,  the  calculation  does  not 

reproduce  the  sharp  decrease  in  stress,  nor  the  subsequent  recovery.  Interestingly, 
the  experimental  record  recovers  to  approximately  the  calculated  values  just  as  the 
gage  location  is  covered  by  the  moving  projectile.  At  gage  #2,  the  experimental 
record  begins  with  a  sharp  decrease  but  then  suddenly  levels  out  to  a  curve  which 
is  nearly  parallel  to  the  computed  response.  The  brief  comparison  at  gage  #3 
ends  when  the  computation  predicts  the  compaction  wave  has  reached  the 

downstream  boundary.  The  transducer  records  beyond  approximately  260  us  may 
have  been  influenced  by  wave  motion  reflected  from  this  downstream  boundary. 

The  computation  suggests  that  a  smooth  expansion  wave  system  propagates 
into  the  compressed  aggregate  as  the  projectile  velocity  slows  down  along  the 
quadratic  path  in  time.  This  expansion  wave  system  provides  an  explanation  for 
the  decreasing  initial  maximum  value  reported  by  successive  transducers,  and  the 
gentle  decay  seen  by  gage  #2.  However,  it  offers  no  explanation  for  the  sharp 
decrease  and  subsequent  recovery  seen  by  gage  #1,  or  the  initial  sharp  decrease 
seen  by  gage  #2.  Sandusky  [53]  has  speculated  that  these  events  may  be  the 
result  of  a  slightly-recessed  4  mm  diameter  transducer  responding  to  a  sliding 
aggregate  composed  of  50  um  diameter  particles. 

C.  Influence  of  Energy  Release  /  PDC  Experiment 

Our  previous  theory  [34]  for  a  quasi-steady  reactive  shock  wave  suggested 
that  ignition  of  the  compacted  aggregate  might  be  controlled  by  heat  transfer  from 
the  hot  gaseous  products  trapped  in  the  pores  (hot  pockets)  of  the  aggregate.  It 
was  postulated  that  ignition  would  begin  when  a  sufficient  thermal  wave  was 
established  in  the  cold  solid  material  surrounding  these  pores.  The  crucial  issue  is 
the  initial  temperature  of  the  gases  in  the  pores.  Recall  that  for  TS-3659 
propellant,  Glancy  et  al.  [15]  conducted  three  PDC  runs  in  Lexan  confinement  at 


62 


impact  speeds  up  to  290  m/s  and  failed  to  see  any  evidence  of  light  on  the 

camera  records,  although  the  tube  ruptured.  This  observation  implies  that 

wave-induced  reaction  in  TS-3659  releases  only  a  fraction  of  the  total  energy,  and 
this  fraction  is  less  than  the  value  which  would  produce  gas-phase  temperatures  of 
say  1500°K.  For  the  computation  of  PDC-80  shown  in  Fig.  31a,  the  amount  of 

energy  released  in  the  wave-induced  reaction  was  adjusted  such  that  the 

temperature  of  the  gases  trapped  in  the  pores  of  the  aggregate  was  approximately 
1450°K.  This  simulation  of  PDC-80  included  wall  boundary  friction,  accounted  for 
heat  transfer  to  the  solid  material,  but  no  further  gas-phase  reaction.  None  of 
the  solid  surface  reached  the  473°K  surface  temperature  ignition  criterion,  and 

hence  no  ignition  occurred.  Since  the  transducers  clearly  show  that  reaction  did 

begin  before  the  end  of  this  calculation,  it  is  concluded  that  heat  transfer  from  an 
inert  hot  gas  trapped  in  the  pores  cannot  be  the  sole  source  of  ignition.  At  least 
for  TS-3659  propellant,  there  must  be  some  additional  source  of  reaction  which 
leads  to  ignition  of  the  aggregate. 

This  study  explored  the  possibility  that  a  delayed  gas-phase  reaction,  which 
converts  R,  to  R^,  is  the  additional  source  of  energy  which  ignites  the  granular 

solid  material.  Baer  and  Nunziato  [19,20]  have  postulated  a  "compressive"  reaction 
with  a  sudden  increase  in  energy  release  which  accomplishes  the  same  result.  The 
discussion  example  is  based  on  PDC-80,  simulated  with  a  constant  projectile  speed 
of  127  m/s  which  gives  a  good  representation  of  the  projectile  locus  up  to  a  time 
of  200  us  (see  Fig.  31b).  Assuming  a  96  us  delay  (and  the  parameters  listed  in 
Table  I)  produces  the  results  shown  in  Fig.  32a  for  gage  # L  In  a  rough  sense, 

the  runaway  pressure  curves  have  the  correct  shape  but  the  "break  away"  is  more 
abrupt  than  the  experimental  data  would  indicate.  This  behavior  is  typical  of  all 
computations  which  assume  that,  immediately  after  ignition,  combustion  produces 
full  energy  release  (solid  ->  R2).  The  two  curves  in  Fig.  32a  differ  only  in  the 
amount  of  surface  area  within  ^the  aggregate  which  is  involved  in  the  combustion 
process.  The  steeper  curve  (chain-dash)  is  produced  by  combustion  of  the  full 

surface  area  of  the  deformed  grains,  while  the  solid  curve  follows  from  combustion 
of  only  the  exposed  surface  area  of  the  grains.  The  "knee"  in  this  curve  in  the 

(text  continues  on  page  70) 


Table  I  -  Input  Parameters  for  PDC-80 
60.2%  TMD  TS-3659  Propellant 

Gas-Phase  Reaction:  t°gj  =  85  microsecs 

Agr  =  1.3x10  6(sec-cm3)-1 
E  gr  =  7500  K 

A egr  =  900  cal/g 

Combustion  (Solid  ->  R?):  r(cm/s)=.244  P(MPa)  0.857 


Energy  Release: 

(a)  Solid  ->  Rg 

(b)  Solid  ->  R  j 

(c)  Rj  ->  Rg 
Wave-induced  Reaction: 


1104  cal/g 
204  cal/g 
900  cal/g 
(b) 


=  i.8% 
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propellant  does  not  ignite. 

Figure  31  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-80: 
160  m/s  impact  of  Lexan  projectile  on  602%  TMD  TS-3659  ball  propellant 


(ram)  aoirs^siQ 


co  co 

p  — 
U 

O  V 

c 

0) 

S-o 

«2  s 
>->  o 
«02 


«4“ 

3  s 

8  . 
•nr* 
O  CNJ 

a 


u 

<2  T3 
a> 


co 


v 

a. 


'In' 

4-» 

G 

U  CO 

02  S 

g 

fc  8 

£ 

•§  i 
§  E 

to 

£  3 

0.3 

0) 

T3  ^ 
0) 

0) 


w  cj 

a  -a 


S 

D 

co  — 

CO 

<c5  ’■p 
a) 

o  w 

G3 

cO 

8 

#C 

sg 

0*  i  o 

So  ♦* 

OQ  4* 
O  CtC 

i 

s-g  > 

W°3 

(M  O 


c 

a> 

s 

u 

<u 

a 

x 

0) 


65 


Figure  31  -  Comparison  of  model  predictions  and  NSWC  data  [15]  from  PDC-80: 
160  m/s  impact  of  Lexan  projectile  on  60.296  TMD  TS-3659  ball  propellant 
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aggregate.  Chain-dash  curve  based  on  combustion  of  full  surface  area  of  grains. 

Figure  32  -  Comparison  of  model  predictions  and  NSWC  data  [42]  from  wall- 
mounted  pressure  transducer  at  38.1  mm  (gage  #1)  for  PDC-80:  160  m/s  impact  of 
Lexan  projectile  on  60.296  TMD  TS-3659  ball  propellant. 


(edK)  3S9J1S  I^TP^  9-mpqH 


o 

in 

o 

in 

© 

i n 

o 

in 

o 

in 

o 

m 

© 

in 

o 

in 

o  in 

o 

£ 

m 

G\1 

o 

t'- 

in 

CM 

o 

r- 

m 

CM 

o 

r*- 

in 

CM 

CM 

■d* 

CO 

n 

n 

cm 

CM 

CM 

CM 

”4 

wM 

1 

«  c  «r 

>  4/  C 
fc*  g*  •  *«* 

3  •>  «3 

O  >  I- 

bo 

12  «fu- 

o  «  ° 

02  60 
£  B3 

1*4  rtj 

bC  S 
«S  “«j 

s-.  (0 

4)  tfl 

o  g*E 

t".  u  3 
O  co 

II  «- 

eu.S  §. 
x 

tf  a  • 

g’g'S 

S  60 

£  S 

CLw-  O 
O  +3 
e  in 

o  «J  2 
•■5  g-2 

1 05 1 

Sg 

+•  2  c 
«S  o 

g  to  Id 

2-  I 
1*3-0 
0) 

•2°s 

ta  -  « 

■*->  5 
2  .2  -a 

0.-M  ® 

c  “j: 
fi  3  w  • 
O-Q  <&  £ 

w  SO  Oh 

"®  8  A 

•s  -u 

£  fi.'*'  M 

J3  U  A  g 

la.=f 


67 


Figure  32  -  Comparison  of  model  predictions  and  NSWC  data  [42]  from  wall- 
mounted  pressure  transducer  at  38.1  mm  (gage  #1)  for  PDC-80:  160  m/s  impact  of 
Lexan  projectile  on  60296  TMD  TS-3659  ball  propellant 
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range  200  -  300  MPa  is  a  typical  consequence  of  the  slower  energy  deposition 
rate  produced  by  combustion  of  only  the  exposed  surface  area.  This  "knee"  results 
from  the  slight  expansion  wave  caused  by  the  formation  of  the  rapid  pressure 
buildup  region  into  a  forward  propagating  wave  system  and  its  movement  away 
from  the  projectile  face. 

Changing  the  energy  release  mechanism  has  a  substantial  influence  on  the 
rate  of  pressure  build-up.  Setting  Ptf=70  MPa  (combustion  of  exposed  surface 
area  of  the  solid  is  accompanied  by  partial  energy  release  while  P  <  Pt  ) 
produces  the  results  shown  in  Fig.  32b.  This  prediction  for  the  gage  w=l  record 
comes  much  closer  to  simulating  the  gentle  lift-off  behavior  seen  by  the 
experiment.  There  are  certainly  no  claims  of  uniqueness  here,  but  none  of  the 

other  schemes  attempted  by  this  investigation  would  reproduce  this  behavior.  The 
two  predictions  in  Fig.  32b  differ  only  in  the  amount  of  surface  area  involved  in 
the  combustion  process  after  P  >  Pt  ;  in  this  example  surface  area  plays  a 
secondary  role,  but  there  are  t&>  few 1 experimental  data  points  to  decide  which 
assumption  is  better.  The  dramatic  influence  of  changes  in  the  value  of  Pt  is 
shown  by  the  results  in  Fig.  33.  The  dashed  curve  is  produced  by  full  energy 

release  after  ignition,  the  solid  curve  (Ptr=70  MPa)  is  repeated  from  Fig.  32b,  and 
the  chain-dashed  curve  assumes  partial  energy  release  up  to  P  =300  MPa  (quite 

unrealistic).  A  faster  gas-phase  reaction  would  tend  to  steeperP  the  latter  curve. 
If  there  were  ever  any  doubt,  these  results  show  the  important  role  played  by  the 
combustion  mechanism  in  the  early  portion  of  the  transition  to  detonation.  And, 

of  course,  the  mechanism  may  change  with  different  materials.  Further  study  will 
be  required  in  this  important  area. 

The  predictions  for  PDC-80  discussed  above  assume  a  constant  projectile 

speed  of  127  m/s.  Re-computing  these  cases  using  a  curve-fit  (see  solid  line  in 

Fig.  31b)  to  the  experimental  data  for  projectile  path  displayed  in  Fig.  31b  results 
in  only  minor  changes  to  the  predicted  gage  response.  This  is  not  true  for  the 

simulation  of  PDC-81,  however.  Figure  34  shows  a  full  simulation  of  PDC-81  using 
the  "dual-linear"  projectile  locus  (Fig.  24a)  which  matches  the  given  experimental 

data.  The  reaction  parameters  are  the  same  as  in  PDC-80,  except  t^  j  =65  us 
and  Ae  °=850  cal/g.  The  prediction  for  gage  #1  based  on  combustion  of  the 

total  si^face  area  of  the  grains  (chain-dashed  curve)  exhibits  a  steep  rise  which  is 

not  supported  by  the  experimental  data.  The  other  prediction  for  gage  #1  (solid 

line)  assumes  combustion  of  only  the  exposed  surface  area,  and  exhibits  a  "knee" 

near  200  MPa  similar  to  the  experimental  data.  The  prediction  for  gage  #2 
begins  runaway  similar  to  the  data,  but  then  becomes  a  strong  reactive  shock  wave 
much  like  that  seen  in  PDC-82.  The  experimental  data  from  both  gages,  however, 
show  a  slower  growth  rate.  The  fact  that  the  record  from  gage  #1  actually 
crosses  that  from  gage  #2  suggests  that  an  expansion  wave  is  propagating  from 
the  direction  of  the  projectile  face.  The  computation  shown  in  Fig.  34,  however, 

assumes  the  projectile  speed  is  maintained  at  173  m/s.  Comparison  between 
experimental  data  and  prediction  strongly  suggests  that  an  expansion  wave  of 
unknown  strength  was  influencing  the  response  at  both  gage  locations.  Since  the 
true  projectile  locus  is  not  known  in  this  case,  the  computation  cannot  be  used  to 
verify  the  reaction  mechanism.  Stated  another  way,  any  prediction  for  PDC-81 
which  demonstrates  a  match  with  the  experimental  gage  records  is  based  as  much 
on  an  estimated  projectile  locus  as  it  is  on  a  reaction  mechanism.  This  is  not  the 
conclusion  this  investigation  sought,  but  anything  else  is  a  misrepresentation. 
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VII.  SUMMARY  COMMENTS  AND  CONCLUSIONS 


(1)  The  equilibrium  stress  state  (Section  II)  which  is  calibrated  with  data 

from  the  quasi-static  compaction  experiment  appears  to  provide  a  reasonable 

description  of  the  compacted  granular  aggregate.  Comparison  of  mode!  predictions 
(Section  VI)  to  data  from  wall-mounted  pressure  transducers  in  the  PDC 

experiment  strongly  suggests  that  the  gages  are  responding  to  the  radial  component 
of  this  mixture  stress  tensor. 

(2)  The  equilibrium  shock  wave  theory  (Section  IV)  predicts  that  the 

fundamental  properties  of  a  compaction  wave  (propagation  speed,  and  downstream 
values  of  stress,  density,  and  porosity)  are  functions  of  the  equilibrium  stress  state. 
This  straightforward  computation  provides  an  important  base-line  to  assess  the 

behavior  seen  in  the  PDC  experiment  When  the  experimental  value  of  wave  speed 
exceeds  this  prediction,  some  mechanism  is  "driving"  the  wave. 

(A)  One  possibility  is  wave-induced  reaction  which  converts  a  small 
amount  (order  of  1  per  cent)  of  solid  material  into  gaseous  combustion  products. 
These  gases  become  trapped  in  the  porosity  of  the  compressed  aggregate,  increasing 
its  "stiffness"  to  support  the  higher  wave  speed. 

(1)  Available  evidence  from  several  materials  would  suggest  that  only  a 
portion  of  the  total  heat  of  reaction  is  liberated  in  this  wave-induced  reaction,  and 
hence  the  gaseous  products  are  intermediate  reactants  which  may  be  involved  in  a 
later  gas-phase  reaction. 

(2)  For  several  materials  of  interest,  the  percent  of  solid  participating 
in  the  wave-induced  reaction  appears  to  decrease  with  increasing  wave  strength. 
This  may  imply  that  phase  change  or  other  nonlinear  physics  are  participating  in 
the  wave-induced  reaction. 

(B)  Another  possibility  is  rate-dependent  resistance  to  dynamic 
compaction.  Modeling  this  effect  (Section  III)  by  assuming  that  the  volume  strain 
rate  is  linearly  proportional  to  the  disturbed  force  balance  cannot  explain  the 
observations  of  the  PDC  experiment  More  complicated  constitutive  behavior  (such 
as  dramatically  different  resistance  forces  to  compression  and  expansion)  were  not 
investigated. 

(C)  It  is  also  possible  that  both  (A)  and  (B)  occur  simultaneously. 

(3)  Various  laboratory  experiments  on  SDT  and  DDT  in  granular  energetic 
material  have  observed  abrupt  increases  in  the  speed  of  propagation  of  compaction 
waves  and  the  leading  edge  of  luminosity.  The  DDT  experiment,  in  at  least  one 
case,  has  evidence  of  a  strong  rearward  propagating  wave.  Some  insight  is 
provided  by  a  simple  quasi-steady  model  (Section  V)  which  envisions  an  increase  in 
reactivity  triggered  by  the  collision  of  two  reactive  compaction/shock  waves. 
Simulation  of  experimental  data  shows  that: 

(A)  After  a  time  delay,  the  onset  of  vigorous  combustion  of  material 
adjacent  to  the  piston  face  will  drive  a  second  compressive  wave  system  into  the 
aggregate  which  has  been  formed  by  the  initial  compaction  wave.  The  leading 
edge  of  this  second  wave  system  propagates  at  a  speed  near  the  local  sound 
velocity. 
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(B)  The  analysis  predicts  that  the  abrupt  increases  in  propagation  speed 
observed  experimentally  would  require  the  collision  to  induce  reaction  in  the  range 
of  296  -  1096.  Transient  combustion  (neglected  in  the  model)  behind  the  leading 
edge  of  the  secondary  wave  system  must  play  a  role  in  the  "trigger"  mechanism. 

(C)  The  wave  pattern  produced  by  the  collision  includes  both  a  combined 
compaction/shock  wave  exhibiting  an  abrupt  increase  in  speed,  and  a  rearward 
wave  which  propagates  back  into  material  compressed  by  the  first  two  waves.  A 
simulated  radiograph  taken  shortly  after  the  collision  predicts: 

(1)  The  density  behind  the  combined  shock  wave  may  appear  nearly 
equal  to  the  quiescent  value.  Wave-induced  reaction  generates  enough  gas  pressure 
to  counteract  (resist)  the  normal  occurrence  of  bed  compaction  caused  by  a  large 
value  of  mixture  stress. 

(2)  A  high  density  ("plug")  region  is  trapped  behind  the  slow-moving 

contact  discontinuity. 

(3)  The  spatial  distribution  of  the  high  density  region  will  exhibit  a 
"reverse  step"  as  a  direct  result  of  the  rearward  propagating  shock  wave. 

(D)  The  model  suggests  that  the  appearance  of  a  rearward  propagating 
wave  may  be  a  direct  indicator  that  a  "trigger"  mechanism  is  responsible  for  the 
abrupt  increase  in  wave  speed.  The  lack  of  any  rearward  propagating  wave  may 
signify  that  the  transition  mechanism  is  controlled  by  nearly  complete  combustion 
which  "drives"  from  behind  until  it  catches  the  leading  compaction  front,  as 
suggested  fifteen  years  ago  by  Bernecker  &  Price  [4)  on  the  basis  of  their  DDT 
experiments. 

(4)  The  transient  reactive  shock  wave  model  (Section  VI)  eliminates  many 

restrictive  assumptions  adopted  in  the  previous  models,  and  accounts  for  transient 
combustion  in  both  gas  and  solid  phases.  Application  to  the  PDC  experiment 

highlighted  three  important  areas: 

(A)  Projectile  Path.  Some  uncertainty  associated  with  the  experimental 

projectile  face  locus  after  impact  was  shown  to  have  an  important  influence  on  the 
predicted  response  of  the  wall-mounted  pressure  transducers. 

(1)  The  "linear  decay"  region  in  the  gage  records  can  be  simulated 

with  a  slight  curvature  adjustment  to  the  projectile  face  locus.  Thus,  deformation 
of  the  Lexan  projectile  face  after  impact  may  have  played  a  larger  role  than 

originally  thought 

(2)  The  lack  of  a  precise  projectile  face  boundary  condition  is  a 
serious  impediment  to  verification  of  possible  reaction  mechanisms. 

(B)  Simulation  of  the  PDC  Experiment  Based  on  Inert  Material. 

(1)  Comparison  of  data  from  recent  PDC  experiments  with  predictions 
from  the  equilibrium  shock  wave  theory  [341  suggests  that  granular  melamine  may 
have  aged  and/or  absorbed  moisture  during  a  ten  year  storage,  and  hence  recent 
compaction  experiments  have  seen  a  "softer"  material. 
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(2)  In  general,  the  model  can  closely  simulate  the  PDC  experiment 

when  the  model  boundary  condition  is  taken  from  the  experimental  locus  of  the 
(non-deforming)  aluminum  projectile.  However,  comparison  of  predictions  for  the 
response  of  the  wall-mounted  pressure  transducers  with  the  NSWC  data  suggests 

that  a  sliding  aggregate  composed  of  small-grain  granular  material  can  sometimes 
lead  to  anomolous  readings  from  a  wall-mounted  transducer. 

(C)  Energy  Release  Mechanisms.  Heat  transfer  from  an  inert  hot  gas 

trapped  in  the  pores  cannot  be  the  sole  source  of  ignition.  At  least  for  TS-3659 

propellant,  there  must  be  some  additional  source  of  reaction  which  leads  to  ignition 
of  the  aggregate 

(1)  The  model  simulations  show  that  a  delayed  gas-phase  reaction 

which  converts  the  trapped  reactive  intermediate  species  into  final  combustion 

products  could  be  the  source  of  ignition  of  the  aggregate.  However,  other  reaction 
sources  cannot  be  ruled  out 

(2)  Comparison  of  model  predictions  with  pressure  transducer  data 

from  the  PDC  experiments  suggests  that  combustion  of  the  solid  begins  with  an 
incomplete  reaction  in  the  flame  zone  (as  in  Fig.  21)  and  then  proceeds  to  complete 
combustion  (full  energy  release)  as  the  pressure  field  rises  rapidly.  The  simulations 
also  suggest  that  combustion  of  the  exposed  surface  area  of  the  compacted 

granular  aggregate  is  controlling  the  pressure  build-up.  Further  study  will  be 
required  in  this  important  area. 

(5)  Unresolved  issues: 

(A)  The  wave-induced  reaction  mechanism  was  not  identified. 

o 

(B)  Sandusky's  [14]  observation  that  ”r  At  =  constant"  was  not  predicted. 

(C)  The  reaction  mechanism  which  controls  ignition  of  the  solid  aggregate 
was  not  determined  uniquely  or  calibrated. 
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LIST  OF  SYMBOLS 


a  =  wave  propagation  speed  in  gas  phase  {see  definition  above  Eq.  24} 

o 

ajg  =  longitudinal  wave  speed  in  solid  phase  {see  definition  above  Eq.  24} 

a^  =  distance  between  deformed  particle  centers  {see  Eq.  14} 

a^  =  intercept  value  in  linear  Hugoniot  representing  TMD  material  {see  Eq.  1} 

A  =  pre-exponential  factor  in  gas-phase  reaction  production  term 
gr  {see  Eq.  16} 

A.  =  modified  longitudinal  wave  speed  in  solid  phase 
s  {see  definition  above  Eq.  24} 

A^  =  slope  of  linear  Hugoniot  representing  TMD  Material  {see  Eq.  1} 

b  =  co-volume  {see  Eq.  10} 

B9  =  coefficient  in  0  (e  A  the  strength  of  the  granular  aggregate  {see  Eq.  9} 
cy  =  specific  heat  at  constant  volume 

C-  =  coefficient  used  in  characteristic  theory  {defined  above  Eq.  24} 

C-  =  coefficient  used  in  characteristic  theory  {defined  above  Eq.  24} 

C~  =  coefficient  used  in  characteristic  theory  {defined  above  Eq.  24} 

dp  =  particle  diameter  {used  in  Eqs.  12  and  13} 

e  =  specific  internal  energy  of  the  gas  phase  {  /gas  volume  } 

6 

e  °  =  reference  value  of  specific  internal  energy  of  the  gas  phase 
g  {  /gas  volume  } 

Ae  0  =  relative  value  {above  e  °}  of  specific  internal  energy 


- -  -  v - -  rr  >  —  I - 

of  reactive  intermediate  species,  R^. 


e  =  specific  internal  energy  of  the  solid  phase  {  /solid  volume  } 

Egr  =  activation  temperature  in  gas-phase  production  rate  term  {see  Eq.  16} 
f(Rp/ap)  =  function  defined  after  Eq.  14 

f ,  =  interphase  drag  force  per  unit  total  volume  [positive  when  u  >  u  ] 

0  {see  Eq.  12}  g 

f  =  wall  friction  force  per  unit  total  volume 
w 

F  =  term  defined  after  Eq.  25 
g 

F  =  term  defined  after  Eq.  25 
s 


81 


g(R  /a  )  =  function  defined  after  Eq.  14 

r  r 

h  =  gas-phase  convective  heat  transfer  coefficient  {defined  in  Eq.  13} 

O 

H  =  collection  of  non-homogeneous  terms,  primarily  from  the  balance  of 
8  gas-phase  momentum  {appears  only  in  Eq.  25} 

Hug  =  collection  of  non-homogeneous  terms,  primarily  from  the  balance  of 
4  solid-phase  momentum  {appears  only  in  Ik}.  25} 

H  =  collection  of  non-homogeneous  terms,  primarily  from  the  derivative  of 
f'  the  equilibrium  stress  state  {appears  only  in  Eq.  25} 

H0g  =  collection  of  non-homogeneous  terms,  primarily  from  both  energy 
equations  {appears  only  in  Eq.  25} 

k  =  gas-phase  thermal  conductivity 

o 

m  =  mass  generation  rate  {due  to  pyrolysis  or  combustion  of  solid}  per  unit 

total  volume  {p  x  regression  rate  x  surface  area;  positive  when  solid  ->  gas} 
s 

m  =  mass  removal  rate  of  gas  species  R.  per  unit  total  volume 
gr  {defined  in  Eq.  16}  1 

M  =  molecular  weight  {sec  Eq.  11} 

n  =  number  of  contact  points  in  compacted  lattice  {used  in  Eq.  14} 

Pj  =  coefficient  in  /3g(es),  the  strength  of  the  granular  aggregate  {see  Eq.  9} 

Pg  =  coefficient  in  /3g(es),  the  strength  of  the  granular  aggregate  {see  Eq.  9} 

Pg  =  Pg(Pg,  eg)  =  static  pressure  in  gas  phase 
Pm  =  mixture  pressure  {given  in  Eq.  7} 

Pr  =  Prandtl  number  in  gas-phase 

Pg  =  Pg(Ps»  ®s)  =  spherical  stress  {pressure}  in  solid  phase 

Ptr  =  transition  pressure:  when  P  <  Ptr,  combustion  of  solid  produces  Rj, 

when*  Pg  >rPtr,  combustion  of  solid  produces  Rg 

q  =  heat  transfer  rate  {from  gas  to  solid}  per  unit  total  volume 

R  =  defined  in  Eq.  2  {related  to  volume  strain} 

Rq  =  universal  gas  constant  {used  in  Eq.  10} 

Rj  =  denotes  reactive  intermediate  gas-phase  species  {see  Fig.  6} 

Rg  =  denotes  final  gas-phase  combustion  product  species  {see  Fig.  6} 

Rp  =  radius  of  deformed  particle  {see  Eq.  14} 
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Re,  Re  =  Reynolds  numbers  defined  in  Eq.  12 

S.  =  exposed  surface  area  of  the  deformed  lattice  per  unit  total  volume 
1  {see  Eq.  14} 

tjjgj  =  "delay"  time  before  onset  of  gas-phase  reaction  {computed  in  Eq.  15} 

T  =  static  temperature  of  the  gas-phase  mixture 

O 

u  =  "particle"  velocity  in  linear  Hugoniot  representation  {see  Eq.  1} 
u  =  velocity  of  the  gas  phase 

o 

ug  =  velocity  of  the  solid  phase 

U  =  u  -  V  =  gas-phase  velocity  relative  to  the  moving  coordinate  system 
8  8  c 

U  =  u  -  V  =  solid-phase  velocity  relative  to  the  moving  coordinate  system 

o  S  v 

VC  =  local  velocity  of  the  moving  coordinate  system 

Vp  =  projectile  velocity  {usually  in  reference  to  the  NSWC  PDC  experiment} 

Vg  =  shock  wave  velocity 

W  =  dz/dt  -  U  {used  in  formulation  of  characteristic  equation  system} 

o  O 

Wg  =  dz/dt  -  Ug  {used  in  formulation  of  characteristic  equation  system} 

Yj  =  mass  fraction  of  reactive  intermediate  species,  Rj 
Yg  =  mass  fraction  of  final  combustion  product  species,  Rg 

a  =  coefficient  used  in  characteristic  theory  {defined  above  Eq.  24} 

als  =  coeff‘cient  used  characteristic  theory  {defined  above  Eq.  24} 

a  a  =  coefficient  used  in  characteristic  theory  {defined  above  Eq.  24} 

0  =  PJe)  =  0  -  P_  =  strength  of  compressed  aggregate  {see  Eqs.  8  and  9} 

S  So  S  g 

7  =  ratio  of  specific  heats  {gas  phase} 

rQ  =  Mie  Gruneisen  coefficient  {assumed  constant  here} 

€  —  solid  volume  fraction  =  solid  volume  /  total  volume 

€  =l-e  =  gas  porosity  =  gas  volume  /  total  volume 
8  8 

5  =  defined  {and  used  only}  in  Eq.  9 

rj  =  fraction  of  incoming  solid-phase  mass  flux  converted  to  gas  phase 

{accounts  for  wave-induced  reaction  in  reactive  shock  wave  solution} 
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Uc  =  compaction  viscosity  {see  Eq.  10} 

|  =  |(ej  =  ratio  of  0  to  a  {see  Eq.  4} 
s  r  8 

p  =  gas-phase  density  =  mass  of  gas  /  gas  volume 
8  =  f.  p  =  mass  of  gas  /  total  volume 

©  5  o 

=  e  p„  +  e„  p„  =  mixture  density 
m  s  s  8  8 

p  =  solid-phase  density  =  mass  of  solid  /  solid  volume 
s 

p  =  e  p  =  mass  of  solid  /  total  volume 

0mr  =  radial  component  of  mixture  stress  {see  Eq.  3} 

0  =  axial  component  of  mixture  stress  {see  Eq.  3} 

or  =  radial  stress  component  in  solid  phase  {positive  in  compression} 

0  =  axial  stress  component  in  solid  phase  {positive  in  compression} 

s 

r*  =  coefficient  in  ($(e  X  the  strength  of  the  granular  aggregate  {see  Eq.  9} 

Jl  S  8 

f(e  )  =  defined  in  Eq.  6 
s 

abbreviations: 

DDT  =  deflagration-to-detonation  transition 
HMX  =  cyclotetramethylenetetranitramine 
LANL  =  Los  Alamos  National  Laboratory 

NSWC  =  Naval  Surface  Warfare  Center  /  White  Oak  Laboratory 
PDC  =  piston-driven-compaction  experiment  (NSWC) 

SDT  =  shock-to-detonation  transition 

SNLA  =  Sandia  National  Laboratory  /  Albuquerque  NM 

TMD  =  theoretical  maximum  density 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3138 _ Date  of  Report  AUGUST  1990 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


s.  fSetsend  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.)  . . 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  Suite,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD  Organization 

ADDRESS _ _ 

Address 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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